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SUMMARY 


Recent satellites have carried long extendable appendages of the 
type classified as thin-walled cylinders of open section for such di- 
verse applications as dipole antennas, gravity- gradient stabilization 
booms, and electric-field probes. Actual flight data show that several 
satellites which deploy these long appendages (or antennas, as they will 
be referred to) exhibit varying degrees of anomalous spacecraft body 
motion while exposed to direct sunlight. NASA/GSFC hypothesis is 
that in many cases the observed oscillations are thermally induced. 

This document presents the equations which describe the coupled 
non-planar transverse and torsional vibrations of such an antenna when 
it is exposed to a directional solar thermal field (simlight). To obtain 
a numerical solution, the equations are specialized to the case of a 
clamped-free cylinder of open section having the nominal character- 
istics of the antennas in orbit. Comparison of results with actual flight 
data shows that this hypothesis leads to conclusions consistent with the 
observed phenomena. The numerical examples show that the problem 
of thermally- induced vibration can be eliminated by literally zipper - 
ing the open section of the antenna, in various fashions, so as to effec- 
tively increase its torsional rigidity. Flight data from the Radio 
Astronomy Explorer satellite (which has four 750-foot stable zippered 
antennas) and the Orbiting Geophysical Observatory, OGO VI, (which 
has two 30-foot stable torsionally rigid electron probes) support this 
conclusion. 
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COUPLED THERMALLY-INDUCED TRANSVERSE PLUS 
TORSIONAL VIBRATIONS OF A THIN-WALLED 
CYLINDER OF OPEN SECTION 


by 

Harold P. Frisch 
Goddard Space Flight Center 


INTRODUCTION 

Prior to the flight of RAE* ** all satel- 
lites having long extendable appendages 
used the basic STEM^ type boom, shown 
in Figure 1, for storage and deployment. 

It is made from a long, thin prestressed 
tape which is rolled on a drum for stor- 
age and when extended takes on a tubular 
shape. Such booms have been used on 
OGOl IV and V as electric field probes, 
on ATS-A,§ OVl-10, GGSE " HI, V and VI 
and several other Navy and Air Force 
satellites as gravity -gradient booms. 

Anomolous Motion of Booms 

On each of the aforementioned satellites, with the exception of RAE, varying degrees of anom- 
alous spacecraft body motion has been observed. This motion is known to be excited only when the 
satellite is in full sunlight and its magnitude has at times been large enough to impair severely the 
performance and stability of the satellite. =•'* 


UNFURLING ELEMENT 



Figure 1 — The STEM principle. 


*RAE — Radio Astronomy Explorer 
tSTEM— storable tubular extendable member 
tOGO— Orbiting Geophysical Observatory 
§ATS— Applications Technology Satellite 
liGGSE— Gravity -Gradient Stabilization Experiment 

**General Electric Study Report for NASA/GSFC Document No. 695D4215, 15 Oct. 1968 '^Gravity Gradient Rod Thermal Flutter Study 
Report” 
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The problem first manifested itself during the initial phases of the flight of OGO IV. Soon 
after the deployment of the 60’ open-section cylinder, used as an electric-field probe on this sat- 
ellite, in full sunlight, a disturbance was noted by the control- system sensors. They indicated that 
the satellite body was oscillating ±1.5° about the axis of maximum moment of inertia at the fundamen- 
tal bending frequency of the satellite boom system. Furthermore, the oscillation persisted whenever 
the satellite was in sunlight and rapidly damped down every time it entered eclipse. From a knowl- 
edge of the elastic modes of vibration of OGO IV it has been calculated* that a ±1.5° oscillation of 
the satellite body corresponded to a ±20-foot oscillation of the antenna tip. 

Before the launch of OGO IV, a study had been completed (Reference 1) which showed that, 
because of the extreme weakness of the boom in torsion, the static -thermal- equilibrium shape 
that the boom could assume in a solar thermal field was not unique. It was postulated^ that, 
if torsional motion were excited, the thermal equilibrium position toward which the antenna 
was moving would be continually changing in time. It was further hypothesized that, if the 
proper relationship existed between the motion of the position of thermal equilibrium and the mo- 
tion of the boom, the observed anomalous motion could be thermally-induced. 

This paper is intended to be a complete documentation of the theory of thermally-induced vibra- 
tion of spacecraft booms; the basic findings have been previously published in a highly condensed form.^^ 

Soon after the initial postulation that the anomalous motion on OGO IV could be caused by a 
thermally-induced oscillation of the attached boom, two other variations of the theory were postu- 
lated by other personnel" of Goddard Space Flight Center. 

The first variation, set down by G. Banks,** showed that a cantilevered beam could be driven 
unstable by solar energy if the source were oriented behind the clamped root. The solution was 
obtained by investigating the stability characteristics of the equation which defined the motion of 
a single-degree-of-freedom analogy of the system. Results of this analysis indicated also that 
only a small amoimt of damping is needed to stabilize the system. This variation was therefore 
dismissed as not practically applicable to present spacecraft problems. 

The second variation, set down by J. Donohue,^^ uses a unique representation of the coupling 
between torsional and transverse displacement to derive a single-degree-of-freedom analogy of 
the system. While some of the assumptions made must be termed gross simplifications of the 
problem, the treatment is a perfect example of how sound engineering judgment can be incorporated 


*GSFC Memo 732-46, H. P. Frisch to H. C- Hoffman, Aug. 18, 1967, "Transfer Function, Tip Deflection and Root Bending Stress of 
OGO-D Boom As Derived from Flexible Body Analysis." 

tGSFC Memo 732-56, H. P. Frisch to H. C. Hoffman, Aug. 29, 1967, "Thermally Induced Vibrations of the OGO-D Boom With Inference 
to ATS.” 

Jj. H. Donohue and H. P. Frisch, "Thermoelastic Instability of Open Section Booms,” Symposium on Gravity Gradient Attitude Con- 
trol, Dec. 3-5, 1968, Los Angeles, Calif. Also, NASA Technical Note TN D-5310, Dec. 1969- 

§H. P. Frisch, "Thermally Induced Vibrations of Long Thin Walled Cylinders of Open Section,” ASME/AIAA 10th Structures, Structural 
Dynamics and Materials Conference, April 14-16, 1969, New Orleans, La. 

II J. Donohue, H. Hoffman, H. Price, B, Zimmerman, and G. Banks. 

•*GSFC Memo, Sept. 1, 1967, G. F. Banks to R. J. Darcey, “Boom Flutter Due to Thermal Lag.” 

ttcSFC Memo 732-194, J. H. Donohue to H. C. Hoffman, "Preliminary Review of Lumped Parameter Simulation of OGO Thermally In- 
duced Boom Oscillations,” Feb. 29, 1968. 
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into a simple model to yield useful and reasonably accurate answers to an extremely complex 
problem. 

Unknown to GSFC at the time, G. Augusti at the University of N^les simultaneously initiated 
an investigation into the stability of elastic structures loaded by non- conservative forces, namely, 
radiant heat.* The problem solved by Augusti was virtually the same as that solved by Banks. 

Both investigators arrived at the same basic conclusions. 

Subsequent to the initial work done by NASA/ GSFC, Y. Y. Yu^ made use of G. Banks* original 
idea and extended it to the continuous-beam case. His results, however, are 180° out of phase with 
the results obtained from the single-degree-of-freedom analogies presented by Banks and Augusti. 
The present author believes that the results of the single-degree-of-freedom analogy are correct 
and that the results of the continuous -beam analysis should be closely approximated by its lumped- 
parameter analogy,^ 

The work to be presented herein derives the equations which define the dynamic response 
of a clamped-free thin-walled cylinder in a solar thermal field. The equations are specialized 
to the case of a cylinder of open section; however, modifications can be introduced which en- 
able one to study cylinders of closed section. Cylinders of zippered section, that is, cylinders of 
open section which have had their seams literally zippered, also can be studied by the proper adjust- 
ment of the torsional stiffness parameters and by the introduction of a reasonable thermal model. 

While it is recognized that the simplifying assumptions to be introduced will prevent describ- 
ing the dynamic response of the boom exactly, it must be emphasized that the purpose of this 
analysis is to define the basic cause of the observed oscillations and to obtain an approximate esti- 
mate of its steady state amplitude. Furthermore, considerable scatter in experimentally determined 
values of stiffness, damping, and thermal characteristics imply that there is no point in an overly 
refined analysis. 

The results obtained by numerical solution of the derived equations show that the postulated 
theory^ (if accepted) leads to conclusions which are consistent with observed phenomena. 

The question of the validity of the theory of thermally-induced vibration as proposed by God- 
dard Space Flight Center to explain this phenomena (and presented herein) was the subject of a 
session at the Symposium on Gravity Gradient Attitude Control sponsored by the Air Force and 
Aerospace Corporation in December 1968. It was the general conclusion at the symposium that 
this theory is valid and does indeed define the cause of the anomalous motion observed on OGO IV 
and other satellites deploying long flexible cylinders of open section. 


•G. Augusti, "Instability of Struts Subject to Radiant Heat,” 12th International Congress of Applied Mechanics, Stanford University, 

Aug. .26-31, 1968; in Meccanica, No. 3, 1968. 

ty. Y. Yu, "Thermally Induced Vibration and Flutter of a Flexible Boom,” AIAA 7th Aerospace Science Meeting, New York, Jan. 20-22, 1969. 
Also, Journal of Spacecraft and Rockets, Aug. 1969* 

^See forthcoming technical comments by Jordan and Augusti in Journal of Spacecraft and Rockets, 

^GSFC Memo 732-56, H. P. Frisch to H. C. Hoffman, Aug. 29, 1967, "Thermally Induced Vibrations of the OGO-D Beam with Inference to ATS." 
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Mathematical Model and Outline of Solution 


The model which is assumed to adequately describe the boom is a long, uniform, continuous, 
thin-walled cylinder of open cross-section. It is assumed to be clamped at one end both in two- 
axis bending and in torsion. The other end is free to warp and deflect and has attached to it a tip 
weight which has both translational and rotational inertia. The effect of the shear center axis being 
displaced from the centroidal axis is assumed to be negligible, and the bending stiffness about any 
axis passing through the centroid of the cross-section is assumed to be constant. 

The system is assumed to be free of all disturbances except those arising from the transla- 
tional and torsional motion of the cylinder relative to an inertially fixed coordinate system, those 
arising from viscoelastic damping both in two-axis transverse bending and in torsion, and those 
arising from the continual attempts of the cylinder to relieve its state of thermal stress by bending 
and twisting toward a position of thermal equilibrium. 

The response of the cylinder to a time- varying non-uniform thermal stress distribution is as- 
certained by applying the techniques of elementary strength of materials. This approach is justi- 
fied by Boley in Reference 2 (Chapter 10). 

Accordingly, the Timoshenko theory of bending and torsion of thin-walled members of open 
section as expounded in References 3 and 4 is directly applied. That is, for small-angle bending 
and torsion the equations are separable and can be solved independently. Furthermore, the stress 
couples that produce pure bending are proportional to curvature, and the stress couple that pro- 
duces pure torsion is countered by the resistance of the cylinder to both torsion and warping. 

The equations of vibration are derived by making use of two coordinate systems; one is fixed 
at the clamped end and the other slides with constant rate along the boom. The use of vector tech- 
niques as discussed by Landau and Lifshitz (Reference 5) and A. E. H. Love (Reference 6) makes 
it possible to derive these equations relative to the sliding coordinate system. 

When long rods are to be studied, special consideration must be given to the development of 
the equations that define the forces associated with the kinetic reactions of accelerating mass ele- 
ments. Since the antennas under study may be of arbitrary length, and non-planar transverse mo- 
tion is assumed to occur, the deflected shape of the boom at each instant of time must be considered 
when such forces are computed. From the computation it is seen that significant coupling exists be- 
tween translational and torsional motion. This coupling defines, for relatively long booms, the 
principal mechanism by which transverse and torsional motion are coupled together. For booms 
of short length, however, this effect becomes of decreasing importance compared to such neglected 
effects as the shear center being offset from the centroid and the non-symmetry of the cross-section. 

With torsional motion existing in a solar thermal field, there will be a continuous change in the 
flow of heat around the perimeter of the boom. The rate at which heat will flow from the illuminated 
side to the dark side is governed by the thermal time constants of the boom. The overall temperature 
distribution and hence the thermal loading is therefore a non-linear function of the boom’s torsional 
response and its thermal time constants. 
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It is assumed that there is zero heat flow along the length of the boom and that the heat flow 
aroimd the perimeter is describable by the standard equation for heat conduction and radiation in 
a one-dimensional solid. Heat is assumed to be radiated uniformly around the entire perimeter, 
and to be absorbed only on the sunlit side according to a cosine law. By solving the thermodynamic 
equation at a sufficient number of thermal stations along the length of the boom, an accurate de- 
scription of the overall temperature distribution, and hence the thermal loading, is obtainable. 

The three equations that define the motion of the cylinder in the transverse and torsional di- 
rections and the thermodynamic equation that is defined at each thermal station are partial differ- 
ential equations of the type which may have their solutions expressed in terms of orthogonal func- 
tions. That is, the solutions to the vibration equations are expressed in terms of their normal 
modes of vibration and the corresponding generalized displacement coordinates; similarly, the 
solutions to the thermodynamic equations are expressed in terms of the normalized thermal modes 
and the corresponding generalized thermal coordinates. 

The actual solutions to be presented are obtained by numerically solving simultaneously the 
equations which define the generalized displacement and thermal coordinates associated with the 
most significant vibration and thermal modes. The solutions show that the theoretically predicted 
amplitudes of the antennas on OGO IV and OGO V are consistent with those observed from actual 
flight data. 

Solutions are included that show how variations in various significant and controllable param- 
eters can affect the steady state amplitude of the thermally-induced vibrations. From these solu- 
tions it becomes evident that the problem of thermally-induced vibrations can be eliminated by in- 
creasing the torsional rigidity of the boom by two to three orders of magnitude. This improvement 
can and has been made by various manufacturers by zippering the seams of the open-section 
cylinder. 

Both RAE, launched September 1968, and OGO- VI, launched June 1969, have torsionally rigid 
booms. As predicted, the four 750-foot booms on RAE and the two 30-foot booms on OGO-VI are 
stable. 

THE THERMODYNAMIC MODEL 

The accurate prediction of whether or not thermally-induced vibrations can be excited depends 
upon the thermal model used to describe the thermal loading. The model to be used must be capable 
of taking into account large-angle torsional motion of the boom relative to a fixed sun line. It must 
also be capable of taking into account the thermal properties of the boom material, which cause 
heat to flow at a finite rate. 

The development of the thermal model to be used will be based upon the fundamental hypothesis 
governing the conduction of heat in solids as discussed by Carslaw and Jaeger in Reference 7. 
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Since the amount of twist per unit length is restricted to small values by the assumptions of 
linear elasticity, and since the temperature gradients across and along the boom are small compared 
to its mean temperature, the assumption that zero heat is conducted along the length of the boom 
can be made and the equation linearized. It is thus possible to determine the entire temperature 
profile by solving the one-dimensional equation of heat conduction and radiation at a sufficient num- 
ber of thermal stations along the boom length. The thermal profile between each thermal station 
is easily determined by linear or higher-order interpolation techniques. 

The question to be answered before this method is used to describe the thermal loading is 
whether or not the detail it provides is necessary. To answer this the thermal response of a unit 
boom element is studied when it is made to oscillate at predefined rates and amplitudes in a solar 
thermal field. The results show clearly that it would be very difficult to construct and justify a 
reasonably simple model approximating the complete thermal profile for arbitrary torsional motion. 

Accordingly, the equation describing the temperature change around a unit element of boom 
length is derived for arbitrary rotational motion of the sun relative to the cross-sectional principal 
axes of inertia. Once this equation has been obtained, it is a simple matter to: 

1. Study the thermal response due to a predefined sun motion. 

2. Study the thermal response at any thermal station along the boom. 

By defining a thermodynamic equation at each thermal station and by solving each simultane- 
ously with the vibration equations so that the sun's relative orientation can be continually updated, 
it is possible to obtain a continuous and accurate record of the thermal loading. This information 
can in turn be used in the vibration equations to update the effective forcing function due to thermal 
loading. 


Thermodynamic Equation of Heat Conduction and Radiation 

Figure 2 depicts the type of thin-walled cylinder of open section which is studied in this re- 
port. It has arbitrary length and radius, its wall thickness must be definable as thin (i.e., it can- 
not support a thermal gradient), and the angular amount of overlap should be not less than 0° nor 
greater than 180°. 

The assumptions upon which the thermodynamic equation of heat conduction and radiation is 
based are identical with those set down in Reference 1 and used to study static thermal bending 
plus twist of long thin- walled cylinders of open section. These assumptions are: 

1. Heat is conducted only in a circumferential direction around the entire perimeter. 

2. The entire cross-section loses heat by radiation. 
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3. The heat absorbed at a point on the sunlit side of the boom element is proportional to the 
cosine of the angle between the surface normal and the sun line. 



4. Heat transfer across the interior of the 
cylinder is assumed to be negligible; how- 
ever, heat transfer from the overlapping 
surface to the overlapped surface is ap- 
proximated by a radiation term because no 
smooth flush contact exists between the 
surfaces. 

While slight differences exist between these 
assumptions, especially in the treatment of the 
overlap region, and those used by Florio and Hobbs 
(Reference 8) for their analytic model, both anal- 
yses lead to results that are consistent with ex- 
perimental findings and with each other. 



Figure 2— Thin-walled cylinder of open section 
(a typical boom segment). 


It follows from the fundamental hypothesis for the mathematical theory of conduction of heat 
in solids (Reference 7) that, if there is no heat generated within or radiated away, the differential 
equation that defines the heat flow around a unit length of the boom shown in Figure 2 is 


^ r}2T(s, Z. t) , r9T(s. z, t) _ ^ 

Kj h - /-'Cli - 0 


where the conditions at the boundaries are 


r9T(S, 2. t) 
d s 


f9T( s , z , t ) 
f9 s 


( 1 ) 


(2) 


and the symbols used are defined as: 


K 


T ~ 


h = 


P = 
c = 


s = 
p = 
z = 
t = 


T(s, z, t) = 


thermal conductivity, Btu/ (sec-in. -°F) 
cylinder wall thickness, in. 
weight density of material, Ibs/in.^ 
specific heat of material, Btu/(lb-°F) 

arc length measured from outer to inner seam around the cross-section, in. 
total perimeter of cross-section, in. 

coordinate used to measure arc length along longitudinal axis of boom, in. 

coordinate used to measure time, sec 

absolute temperature at the position s, z at time t, °R. 


If it is assumed that the boom radiates as a black body and that the heat input can be defined 
by the function T(s, z, t) then the thermodynamic equation of heat conduction and radiation is 
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given by 


K^h 


f( s, z, 
ds^ 


t) aT(s, z. t) 

pch 


crel^ (s, Z, t) - T(s, 2 , t) 


(3) 




where 


= Stefan- Boltzmann constant, 

Btu/ (sec-in. 2 

€ = emissivity, 

T(s. z, t) = heat input at position (s, z) and 
at time t, Btu/(sec-in.^), 

The heat input function T(s, z, t) can be 
written as a function of the sun’s orientation 
relative to the cross-sectional principal axes 
of inertia. Using Figure 3 as a guide, we let 

SL(z, t) = unit vector emanating from the 
geocenter of the cross-section 
to the sun’s relative position at 
(z* t) „ 

n(s) = unit vector normal to the surface 
of the boom at position s 


Figure 3 — Boom cross-section at station z. 


-(s, z, t) 


if SL(z, t) ■ n(s) > 0 

. (4) 

if SL( z . t ) ■ nr s ) < 0 


(this function is used to distinguish mathematically between the sunlit and shaded sides of the boom), 


jg = solar radiation intensity, Btu/(sec-in.^), 
= surface absorptivity, 
r = radius of cylinder, in., 


and define the heat input function by the expression 


J 3 0 .^ -(s, z. t) SL(Z, t) • n(s) 


for 0 < s < 2 t 7 r 


T(s, z. t) = < ^ • ( 5 ) 

( s “ 277T , 2 , t) for 277r £ s < P 

This function is analogous to that used in Reference 1 to study static thermal deflection, and 
yields results which are consistent with the experimentally obtained data reported in Reference 8. 
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Making use of the fact that the maximum temperature difference between any two points on 
the boom surface for any sun orientation is small compared to the steady state mean temperature, 
the thermodynamic equation of heat conduction and radiation (Equation 3) can be linearized. Let 


T(s, z, t) - Tq + T(s, z, t) , 


( 6 ) 


where 

Tq = steady state mean temperature of cross-section, °R, 

T(s, z, t) = temperature deviation from 

and express (s, z, t) as a truncated binomial series expansion about Tq. A direct substitution of 

(s, z, t) ^ + 4 Tq3t(s, z, t) 


and Equation 6 into Equation 3 leads to the one-dimensional linearized thermodynamic equation of 
heat conduction and radiation along with its appropriate boundary conditions. That is, 


(9 ^ T(s, z , t ) 


pc t9T( s , z , 


t) 


K^h 


T(s, 


t) 


h 


h 


T(s, 2 , t) 


(V) 


where 


(9T(s, z, t) 
d s 


s= 0 


(9T(s, z, t) 
ds 


5 = P 


0 - 


(8) 


This equation can be solved for arbitrary sun motion by expressing the solution in terms of 
the orthogonal functions which satisfy its homogeneous part and their associated generalized 
coordinates. 

By defining ^(z, t) to be a function of the torsional motion of the boom at thermal station z, 
a time history of the thermal loading at z can be obtained. To determine the thermal loading along 
the entire length of the boom, ^(z, t) must be defined, and the thermodynamic equation must be 
solved, at a sufficient number of thermal stations to obtain a piecewise continuous description of it. 


Solution to Linearized Thermodynamic Equation 

The linearized thermodynamic equation of heat conduction and radiation, 

z, t) pc aT(s, z, t) _ _ o-eTJ* - T(s, z, t) 

at - K^h T(s, z, t) - 


( 7 ) 
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may have its solution expressed in terms of the orthogonal functions that are solutions to its homo- 
geneous part and that satisfy the stated boundary conditions 


^T(s, 2, t) 
(9s 


(9T(s, z, t) 
d s 


( 8 ) 


The orthogonal functions that satisfy the homogeneous equation 


a ^ T( s , z , t ) pc aT(s, 2. t) 


d s " 


dt 


iqir z. t) 


(9) 


may be obtained by separation of variables. If a solution of the form 

T( s , z , t ) = ' T( s ) 


( 10 ) 


is assumed and substituted into Equation 9, a standard-type eigenvalue problem is obtained. That 
is, the function T(s) must be a solution to the equation 


T(s) 
(ds ^ 


- n_2T(s) = 0 


( 11 ) 


and satisfy the boundary conditions 


dT(s) 

<ds 


dT(s) 

ds 


0 , 


(12) 


where 


4aeT ^ 

“t ■ K.j.h ^ 


(13) 


For a non-trivial solution to Equation 11 to exist, it can be shown, by substituting its general 
solution 


n_s -D 

T(s) = Dj e ^ ^ 


(14) 


into the boundary condition equations, that 0^ must be a solution to the characteristic equation 


n._p -n^p 
e ^ ~ e 


( 15 ) 
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and that the integration constants and D 2 must satisfy the relation 


D, - D 2 = 0 . (16) 

It is easily seen that the characteristic equation, (15), is satisfied for all values of given by 

0^ - ' where n = 0, ±1, ±2, *** , (1^) 


while the eigenvector 


T(s) = T^(s) 


r fi s -0 s' 

Di, [e " 


associated with the eigenvalue is a solution to the equation 


ds ^ 




(18) 


(19) 


and satisfies the boundary conditions 


ds 


s = 0 


ds 


s = P 


0 . 


( 20 ) 


The constants of integration are arbitrary, but can be made unique by requiring the eigen- 
vectors to satisfy the orthogonality relationship 


P 


r 


T^(s)T„ (s)ds 


m, n 


where 


s 


m, n 



/ n 


One is thus led to the expression 


Tn(0 - 


_ nrr 
j2 cos s 


for n - 0 

for n=l,2,3, ••• 


(21) 


( 22 ) 


(23) 
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These eigenvectors, referred to as the "thermal modes" of the system, define the set of orthogonal 
functions which will be used to solve the non-homogeneous equation (Equation 7). 

The values of \ that satisfy both Equations 13 and 17 are given by 


pch ^ pc yP y 


(24) 


and are referred to as the "thermal decay constants." The quantities defined by 


r 


n 



(25) 


are referred to as the "thermal time constants" and define the time it takes the nth thermal mode 
to decay 63.21% in the absence of thermal input. 

These functions and relationships developed with respect to the homogeneous equation, (9), 
can be used as a basis for solving the non-homogeneous equation of heat conduction and radiation, 
( 7 ). 


Assume a solution to Equation 7 of the form 


T(s, z, t) - (z, t) (s) , 

[ 1=0 


(26) 


where 

(z, t) - the time-dependent generalized thermal coordinate at thermal station z and time t 
associated with the nth thermal mode (s). 

By directly substituting the assumed solution into the non-homogeneous equation, (7), and by mak- 
ing use of the functional relationships given by Equations 13 and 19, an infinite series representa- 
tion of Equation 7 can be obtained. That is, 

dq, (z, t) ) (s) T(s, z. t ) - 

dt ■ '^n P “ pchP ’ (27) 

n = 0 ^ ^ 



Since the thermal modes (s) are orthogonal functions and are normalized according to Equa- 
tion 21, Equation 27 can be decoupled to obtain a set of linearly independent equations that define 
the generalized thermal coordinates. These equations are derived by multiplying both sides of 
Equation 27 by C s) and integrating around the cross-section. The resulting equations for 
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n 


= 0, 1, 2, • • • are: 


dq„ (z. t) 

dt - Q„(Z, t). 


(28) 


where the "generalized thermal input" (z, t) is defined by 


Qn (z. t) 



t)-creTJ*] T„ (s)ds 
pchP 


(29) 


Generalized Thermal Input 

The generalized thermal input (z, t) associated with the nth generalized thermal coordinate 
( z , t ) was defined in the previous section by the equation 

I [T(s, z. t)-c.eT„‘*] T„ (s)ds 

Qn pchP ' (29) 


The heat input T(s, z, t) was assumed to be adequately described by the equation 


T(s, z, t) 


J 5 H( s , z , t ) SL( z , t ) • n( s ) 

for 

0 < s <2ttv 

oea^ T'* (s ~2Trr, z, t) 

for 

2'nv < s < P 


(5) 


where the sun-line vector SL(z, t) is an arbitrary but definable function of time at any particular 
thermal station z along the boom. 


Figure 4, showing the position of the sun- line vector SL(z, t) relative to the principal cross- 
sectional axes of inertia at the thermal station z, makes use of the following symbols: 


0 = angular amount of overlap 


Y-axis = cross-sectional principal axis of inertia parallel to the axis of symmetry and directed 
away from the overlap area 


X-axis = cross-sectional principal axis of inertia normal to the Y-axis and passing through the 
centroid of the cross-section, directed in the direction shown 

0(z. t) = angle between the negative Y-axis and the sun- line vector at thermal station z and 
time t. 0(z, t) is measured positive in the direction of increasing s. 
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From these definitions and the geometrical 
relationships shown in Figure 4 it can be seen 
that 


SL( z, t)*n(s) - cos 


[^(z. t)- 7 + l] , 


(30) 


and thus Equation (5) may be rewritten as 

t) -7 


T(s, z, t) = < 


4> 

S(s, z, t) cosIsiCz, t) - 7+7 
for 0£s_<277r 

[To'* +4To3T(s-27rr. z, t)] 
f or 277r < s < P 


(31) 


Figure 4—Orientation of boom cross -section 
at station z relative to sun. 


where 


H(s, z, t) - < 


r 

1 if cos 

0 if cos 


s 0 

S^(z, t) - 7 + -2 


> 0 


0(Z, t) 


s 4 > 

7 "T i 0 


By introducing the following definitions, 


B,„(z,t) - J S(s, z, t) cos |^(z. t) - 7 + -fj (s) ds . 


r 

•'0 


( s ) ds 


r 

^2 tt r 


T fs)ds 


(32) 


(33) 


(34) 


(35) 


= r 

nn I 

^277 r 


T„(s - 2 t 7 t)T„ (s)ds 


(36) 
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the generalized thermal input Q„ ( z . t ) may be rewritten as 


Q„(z. t) 



(z. 




To^in +4 2 ^ q„(z, 

_ m= 0 



( 37 ) 


By recalling that 


T(s, z, t) - Tq + Y (z, t)T^(s) , 

m= 0 

one sees that is proportional to the heat that would be radiated by the overlap if the entire 

overlap had a uniform constant temperature equal to the mean Tq, and that 


m= 0 

is a correction factor which accounts for the actual non -uniformity of the temperature distribution. 
Since the mean is two orders of magnitude greater than the maximum temperature deviation 
from Tq, one may write 


ToE., + 4 


z: 


E2mn ^iT 


( Z. t ) 




( 38 ) 


The generalized thermal input is thus independent of all generalized thermal coordinates and can 
be expressed as 


Q„ (z, t) 


pchP 


( 39 ) 


and it will be shown that (z, t ) can be written as an analytic function of ^(z, t ) , 

It is therefore possible to arrive at an analytic expression that defines the quantity (z, t) for 
an arbitrary piecewise continuous function of the angle ^(z, t) . 

It is recalled that t) is the angle between the sun line and the cross-sectional axis of sym- 
metry and hence is directly proportional to the magnitude of the torsional oscillation at thermal sta- 
tion z and time t. 

By defining 0(z, t) to be a particular function of time, the thermal response of a unit element 
of boom length oscillating, relative to a fixed sun, about its longitudinal axis can be determined. By 
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defining >/'(z, t) to be proportional to the torsional oscillation of the vibrating boom, the thermal 
loading at thermal station z can be determined as a function of time. 

hi addition to the above expression which defines the generalized thermal input Q„ (z. t), it 
will be shown that, in order to arrive at the expressions which define the thermal loading com- 
pletely, the derivative of Q„ (z, t) with respect to the sun angle >/'(z, t) must also be determined. 
That is. 


t) _ J, g, aB,„(z. t) 

(9</;(z, t) pchP (90(2, t) (40) 


is defined so that the equation 


t)\ 

dt \ (90(z, t ) / 


+ \ 


n I (9 0( Z , t ) 


t )\ dQ^ (Z, t) 


^0(Z. t) 


(41) 


(to be derived subsequently) can be solved. 

Recalling that the thermal mode shape (s) is given by 


Tn(0 


for n - 0 


^ n n 

y2 cos ^ s for n-1,2,3, 


enables one to evaluate the integral equations defining 


(23) 


Bln (z- t) • , Ej„ , 


^Bln(^- t) 
di/>(z, t) 


as analytic functions of z , t ) . A summary of these expressions follows : 


^20 B 


= 0 n = 1 . 2 . ■ 


E,o = r<j& 


(42) 

(43) 

(44) 


E 


1 n 


I^P 2h 772 r 

" 1m P " “ 1, 2, • • • . 


(45) 
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For n = 0 ; 4 >> 0 ; - cj>/ 2 ± 4 >{z, t ) < 2?? - 4>/2 : 


Bjo (z. t) - 2 r 

( 9 Bjo(z, t) 
d^z, t) - 0 


For n - 2 ', 4 > =0; ~ 4 >/ 2 <\p(,z, i)< 2 n-<p /2 : 


B.J (z, t) 


t /2 ITT 

2 — cos>//(z, t) 


(z. t) 
( 90 (Z, t) 


}/2TfV 

2 sin i//( z , t ) 


For n - 1, 2, * • . ; 0 >0 ; - ^/2 < 0(z, t ) < tt/2 ~ 0/2 ; with the exception of n = 2 for 0 = 0: 


rp2 


|sin j^C 


z, t ) + + cos 


0 7t" 

0( z , t ) + -2 + “2 


r, 0 3 t7‘ 

+ cos -Y~ [0(z, t) + -J + -Y 


1 ■ f/ V 2n?72 

- -sin 0(z, t ) + -^ cos — p— 


n77r f 0l 2n?7 ^ 

p cos 0(z, t)+ sin — p — 


- (n?7r 




0 ( Z , t ) + 


n77T rmr f if 

p sin 0(z, t) + "2 + “2 


n77r n?7 r r 0 3 t 7 

-p- sin -p- 0(z, t )+-2 + -^ 


0 

0( z , t ) + "2 


0l 2n77- 2 
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For n= 1, 2, tj/2 - 4>/2 1 >A( z . t ) < 3 t7/2 - 4>/2 ; with the exception of n = 2 for ^ = 0 


For n = 


t) 


<?B,„ (z. t) 
50(z, t) 


Bln 


2 i2 rP* mrr f , i/>1 r 

p2_'(n^r)'2 ~P“ • 


(z. t) 
<9i//(z, t) 


2 /2n7zr^ P 

P^ - (ri77r)^ 


n77r 
sin p 


r <^1 

['/'(z. t) + -2 


rm^ r 

COS 2p 


L, 2, • • ■ ; <^> > 0; 3n/2 - 4>/2 <'P(z, t)<2n - 0/2; with the exception of n = 2 for 0 = 0 = 


y2 rp2 


p2 -(n77r)^ 


|y-(z, t) + 


+ COS 


rrrTT [ 3tt 

— p- |^0( z , t ) + “2 ^ 


nTTr 

+ cos “p” 




5] 


“ sin 




2n77^ r 

cos — p 


ri77r 

p“ cos 


r 

1 2n77^ rl 

j^0(z, t) + -2 

J sin p 1 


l/2 rp2 


(rr/ir)^ 


f r 

< cos 0( z , t ) + ~2 


n'nr rtrrr [ 4> 3n 

- sin -p- ^(z, + ~2 ~ ~2 


077 r 077 r ttI 

- sin 0(Z, t) + “2 


- cos 


4^ 

kP(z, t ) + j 


cos 


2n77^ r 
P 


077 r 

+ sin 


U( z , t ) • 


s in 


2n77^ 

P 
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Thermal Loading Associated With Non-Uniform Temperature Distribution 


Since all booms of practical interest are made of materials that have finite coefficients of 
thermal expansion, a non-imiform temperature distribution aroimd and along the boom will give 
rise to a non-uniform thermal stress distribution. This condition the boom will attempt to relieve 
by bending and twisting to a position of static thermal equilibrium. I 

As previously stated, the deflection analysis of beams subjected to thermal loading is per- | 
formed on the basis of elementary strength of materials theory. From this basis it is shown in 
References 1 and 2 that the thermal loading may be adequately described by a resultant thermal 
bending moment distribution that produces pure bending and a resultant thermal torque distribution 
that produces pure torsion. 

Both the thermal bending moment and the thermal torque at any thermal station z can be de- 
fined in terms of the generalized thermal coordinates. The corresponding distributions can be 
obtained by numerically differentiating with respect to z . 

It was shown that the solution to the thermodynamic equation of heat conduction and radiation 


h 


(9 ^ f( s, z . t ) 


t?T(s, z, t) ^ 

pch cjeT^ ( s , z, t ) 


T(s, z. t) 


(3) 


was given by Equation 6 combined with 26: 


T(s, z. t) - Tq + 2]^ (z, t) (s) , (46) 

n=0 

where the generalized thermal coordinates (z, t) can be obtained by solving the ordinary differ- 
ential equation 

dq„ (z. t) 

dt - Qn(Z’0- (47) 


The magnitude of was not specified, but Tq was defined to be the steady state mean tempera- 
ture of the cross-section. Its magnitude can therefore be derived from the expression 


lim q (z, t) 

t-OD 


Qo t) 


0 . 


(48) 


Since 


2J. 

Qo '•) pchP 


(49) 
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and 


^ 0 • pch 


(50) 


it follows that 


2rJ, a 

^4 = ? ! 

® ere [p - r0] 


(51) 


In Reference 1 and in Appendix A it is shown that the thermal bending moment components at 
thermal station z about the local cross-sectional principal axes of inertia are given by 


BM^ (z, t) = e 


,Eh f T(s, 

Jo 


t) (s) ds 


(52) 


r 


BMy (z, t) = -e^Ehj T(s,z, t)x^(s)ds 


(53) 


and that the thermal torque about an axis normal to the cross-sectional plane at thermal station z 
is given by 


Tsc ("• t) 


dV(z, t) d\p(z, t) dqp(z, t) 
d<p(z, t) (?<p(z, t) dz 


(54) 


where 

BM^ (z, t) = component of the thermal bending moment vector in the X-direction, i.e. normal 
to the axis of symmetry (see Figure 4), 

BMy (z, t) = component of the thermal bending moment vector in the Y-direction, i,e. parallel 
to the axis of symmetry (see Figure 4), 

= coefficient of thermal expansion, 

E = Young's modulus of elasticity, 

(s) = distance of the point s on the cross-section to the Y-axis; 


(s) 



(55) 
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Yc (s) = distance of the poiat s on the cross-section to the X-axis; 


Yc (s) - “ r cos (7 - 4) * (56) 

8^ = distance between the centroid and the geometrical center of the cross-section; 

. 4 > 

r sin “2 

Sc = W' (57) 

77 + y 


cp(z, t) = angular amount the cross-section at thermal station z has rotated at time t from 
its zero position, 

V(z, t) = 


sin -2 +^J 


T( s , z , t ) ds 


cos ^7 - + ij J T(^, z, t ) d^ds 


(58) 


= distance between the shear center and the geometrical center of the cross-section; 


4 > 

^ 2 


, ^ 1 . , 
77 + '2-'2Sin0 


(59) 


A direct substitution of Equation 46 into Equations 52, 53, and 54 leads to the equations which 
define the bending moment components and the thermal torque as infinite series. That is, 


BM^ ( z > t ) - Eh 


r CO 

I Tj + 2^ (z, t) (s) 


y, ( s) ds , 


(60) 


BMy (z,t) - -e^Eh 


r 

-'0 


n = 0 


( S ) ds 


(61) 


T.c t) 


( s ) ds 


Jr . pi r o 

'cEhr^l^e^s.n 2 +^J J Tn 

^ ^ n= 0 

■ [ ' '] i 12 Tn (O 


ds 


ai/-(z, t) d(p(z, t) . . 

(9tp( z , t ) dz ’ ’ 
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where 



By introducing the definition of the thermal mode shape ( s ) given by Equation 23 and carry- 
ing out the obvious integrations, one is led to the following; 


(z, t) 


2 Ehr2p2 sin 



n= 2 
n eve n 


(64) 


BM^ (2, t) 


“ 2 /2 Ehr^ cos 



n= 1 
n odd 


(65) 


(z, t) = 2l^e^EhrP< 


r 

n77Pre^ cos 

- 





5qn (Z' t) 

^ Z , t ) drn( Z , t ) 


- (nrrr)2 

077 

<?'/'( Z, t) 

^^cp(z, t) dz 

odd 


— 

J 



( 66 ) 


It is thus possible to obtain infinite series expressions for the thermal bending moment com- 
-ponents and the thermal torque in terms of the generalized thermal coordinates which can be readily 
obtained at every thermal station. 

Special note is made of the fact that these quantities are derived relative to the principal cross- 
sectional axes of inertia at thermal station z. 


Thermal Response of Unit Boom Length to Predefined Torsionai Osciiiation 

Before an attempt is made to define the thermal loading along the entire length of the boom and 
include this thermal loading in the vibration equation, a prerequisite study will be made. 

Two fundamental points must be studied. That is, we must first determine the number of 
thermal modes needed to adequately define the thermal loading, and must then show that the thermal 
model as defined does indeed predict a physically realizable representation of it. 

With regard to the first point, it suffices to state that by comparing results obtained using 2, 

4, 6, and 20 thermal modes (excluding the n = 0 mode) it has been found that the use of 4 modes is 
optimum from the standpoints of accuracy and computer running time. 

The second point is studied by investigating the change in the thermal load on a unit element 
of boom length due to a predefined change in the sun position with time. 

The nominal parameters used to define the geometric and physical properties of the boom ele- 
ment are all obtained from the extensive experimental work carried on at NASA/GSFC on booms of 
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Table 1 


Antenna material 

Perimeter of cross-section 

Total overlap angle 

Wall thickness 

Weight density of material 

Youngs s modulus 

Thermal conductivity 

Specific heat 

Solar radiation intensity 

Thermal expansion coefficient 

Emissivity 

Absorptivity 


= beryllium copper, silver-plated 
= P = 2 in, 

= = 90 deg. 

= h = .002 in. 

= P = .2714 Ib/in,^ 

= E = .19 X 10® lb/in.2 

= 4.167 Btu/(hr-in.-°F) 

= c = .1 Btu/(lb-°F) 

= Jg = 3.065 Btu/ (hr-in. ^ ) 

= - .104 X 10'^ in./(in.-°F) 

= e = .035 



flight quality. These parameters are listed in 
Table 1. 

The hypothetical case of a spinning sun is 


studied in Figure 

5. 

The curves labeled A 

and C correspond 

to 

the particular cases 

fined by 



^(2. t) 

= 

.05 t 

0(Z, t) 

- 

. 50 t 

0(z- t) 

= 

5. 0 t ( radians) 


respectively. The change in the thermal load is 
illustrated by plotting the components of the ef- 
fective thermal bending moment vector (BM^ vs 
BMy as defined by Equations 64, 65) after a 
steady state response has been achieved. 



A. = 0.05f 

B. 0.5 

C. 'A (z,t) = 5.0t 


Figure 5 — Variation in thermal bending moment for 
three spinning sun cases. 


As expected, for the slow- spin case the 
thermal model does yield results consistent 
with those obtained in the static analysis. This 
fact may be verified by comparing curve A with its analogous counterpart in Reference 1. As spin 
rate increases, a temperature-averaging or ”barbecuing^’ effect is to be expected and is evidenced 
by curves B and C. The offset point about which the three curves appear to be centered is a result 
of the fact that the overlapped portion of the cylinder is never exposed to direct radiation. 
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Figures 6 through 14 are included to demonstrate the need for a rather detailed thermal model 
in this analysis. 

In the actual problem, the boom oscillates torsionally in the directional thermal field. To study 
the steady state thermal response of a torsionally oscillating unit element, the sun’s relative orien- 
tation is described by 

0(z, t) = v/;* sin (cD* t - (radians) 

for particular values of , cu*, and 0 ^. 

As in Figure 5 , each curve corresponds to the motion of the tip of the effective thermal bend- 
ing moment vector after a steady state thermal response has been achieved. 

On all figures from 6 through 14 , the curves labeled A, B, C, and D correspond to 0 g = 0 , 77/2, 
77, and 377/2 respectively. 

Table 2 Lists the respective magnitudes of and 0* used in each figure. 




BMy 

A. >//( z, t ) = 7 t/ 4 sin ( 0.5 1 - 0) 

B. '^( z,t)= 77/4 sin (0.5f - 7 t/2) 

= 7t/ 4 sin (0.5t -77 ) 

(z,t)= 77/4 sin (0.5t - 3 7 t/2) 

Figure 6— Variation in thermal bending moment for ' Figure 7 — Variation in thermal bending moment for 
oscillating sun cases. oscillating sun cases. 




BMy 


A. 0(z,t)= V4sin {.05t - 0) 

B. 0 ( z, t ) = 77/4 sin (,05t - V2) 

C. 0(z,f)= 77/4 sin {.05t -77) 

D. 0 (z,t)= 77/4 sin {.05f - 37 / 2 ) 
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The values of cj* listed bracket the frequency range in which nearly all booms of interest will 
respond torsionally. It is interesting to note that outside of this frequency range a much simpler 
thermal model could probably be justified. 

To intuitively justify the results shown, one must determine the magnitude of the thermal time 
constants For the parameters listed in Table 1 these values are 

= 330. 38 sec , 

T j = 9 . 236 sec , 

= 2. 359 sec , 

T ^ - 1 . 052 sec , 

T ^ = . 593 sec . 

A few general statements can be made concerning these figures. 

1. The effect of thermal time constants cannot be neglected in any analysis for which vibra- 
tional response periods and thermal time constants are of the same order of magnitude. 

2. In any particular case of a torsionally vibrating cylinder of open section, the amplitude of 
vibration will vary from station to station along the length. By comparing Figures 6, 9, 12, 
Figures 7, 10, 13, and Figures 8, 11, 14, it is evident that the thermal loading can be quite 
different at different places along the length. A sufficient number of thermal stations along 
the length will therefore have to be used to describe the thermal loading accurately. 

Figures 5 through 14 deal exclusively with the effective thermal bending moment components 
BM^ (z, t) and BMy (z, t). Figures 15 through 24 deal exclusively with thermal torque. To do so, we 



Figure 15 —Variation in thermal torque coefficient for Figure 16— Variation in thermal torque coefficient for 
spinning sun cases. oscillating sun cases. 
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Figure 17— Variation in thermal torque coefficient for 
oscillating sun cases. 



Figure 19— Variation in thermal torque coefficient for 
oscillating sun cases. 



A. ( z, t) = V2 sin ( 5t - 0 } 

B. «//(z,r) = V2 sin (5f - V2) 

C. i//(z,f) - V2 sin (Sr-TT) 

D. i//(z,r) = '^/2 sin (5t- 3*^/2) 

Figure 21— Variation in thermal torque coefficient for 
oscillating sun cases. 



A. </;(z,t) = sin (5t ~ 0) 

B. y\f{z,f) = sin (5f -- V2) 

C. i//(z,t) = ^/4 sin (5t -tt) 

D. t//(z,f) = V4 sin (5t - 3rr/2) 


Figure 18— Variation in thermal torque coefficient for 
oscillating sun cases. 



A. i// (z,r) = ^/2 sin (0,5t - 0) 

B. rp (z,t) = V2 sin (0.5t - V2) 

C. (z,t) = ^/2 srn (0.5t - 7 t) 

D. iff (z,t) = V2 sin (0.5t - 37r/2) 


Figure 20— Variation in thermal torque coefficient for 
oscillating sun cases. 



Figure 22 — Variation in thermal torque coefficient for 
oscillating sun cases. 
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Figure 23— Variation in thermal torque coefficient for Figure 24— Variation in thermal torque coefficient for 

oscillating sun cases. oscillating sun cases. 

plot the temperature-dependent part of (z, t). That is, in all figures from 15 through 24, 


dW( z, t ) 

Ehr vs t) 


is shown for the same respective functions of 0(z, t) used in Figures 5 through 14. 

These figures (15 through 24) are included to provide the reader with a graphic showing of how 
thermal torque is dependent upon relative sun orientation. These figures will play an important 
role in the determination of whether or not a thermally- induced torsional instability can be excited 
for a particular sun orientation. 

The same general statements made previously concerning Figures 5 through 14 can be applied 
here in discussing Figures 15 through 24. Again the reader is directed to Reference 1 for additional 
work pertaining to the static analysis of thermal bending and twist. 

Preliminary test data obtained by R. Predmore and C. Staugaitis at NASA/GSFC indicate that 
the thermal model used herein does tend to yield results which are consistent with experimentally 
observed bending and torsion. More testing is necessary, however, before any direct comparison 
between experiment and theory can be made. 


EQUATIONS OF BENDING AND TORSION 

In order for the results of this work to have practical significance, a certain amount of general- 
ity must be reatined in developing them so that they may be applied to cylinders of open, closed, and 
zippered cross-section with only minor adjustments of the torsional stiffness parameters and the 
thermal model. Since the actual booms used in flight are not mathematically perfect, any attempt to 
describe analytically the non-uniformities in such parameters as bending stiffness, cross-sectional 
diameter, shear center offset, initial curvature frictional damping, etc., leads to equations which are 
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and specialized to a particular experimentally studied boom. The inclusion of these details tends 
to mask the basic mechanism producing thermally- induced vibrations behind a veil of second-order 
perturbations. It is argued that, if the simplified model can be shown to be strongly unstable, it 
would be virtually impossible to show that the inclusion of second-order terms could reverse the 
situation and analytically stabilize the system. 

Accordingly, the boom is modeled as a beam clamped in both two-axis bending and torsion at 
its root. The tip is free to translate and rotate, and has attached to it a tip weight which has both 
translational and rotational inertia. The tip weight is fixed to the boom tip in a way which allows 
it to be free to warp. For the long booms of interest, it is assumed that the effect of the shear 
center axis being displaced from the centroidal axis is negligible and that the bending stiffness 
about any axis passing through the centroid of the cross-section is constant. 

The determination of deflection imder thermal loading is performed under the assumption that 
the longitudinal thermal strain in any fiber passing through an arbitrary cross-section along the 
length of the rod is a linear function of the distance from the fiber to the neutral axis of the 
cross-section. 

Furthermore, the stress couples that produce pure bending are proportional to local curvature, 
and the stress couple that produces pure torsion is a fxmction of the resistance of the cylinder to 
both torsion and warping. By var3dng the magnitudes of the torsional stiffness parameters, it is 
possible to show how the response of the boom changes as the torsional stiffness is increased from 
the case of an open-section cylinder to that of a closed-section cylinder. This transition region in 
effect simulates the torsional stiffness characteristics of zippered-section cylinders. 

As indicated in the Introduction, the small-angle assumptions made are valid over any short 
segment of the boom; however, when extended over very long lengths they do not completely define 
the situation. It is thus desirable to develop the equations of elastic equilibrium from elementary 
considerations. 

This is done by using vector techniques to derive the equations of elastic equilibrium relative 
to a coordinate system which slides along the longitudinal axis of the boom with uniform constant 
rate, relative to an inertially fixed system. 

The body forces associated with the translational and rotational acceleration of mass elements 
along the boom, the dissipative forces associated with viscoelastic damping, and the thermal load- 
ing induced by the thermal environment can also be expressed in vector form relative to either the 
sliding or the fixed coordinate system. 

These equations can, by means of the transformation matrix, be expressed in a common co- 
ordinate system and reduced to a manageable form by the deletion of terms that are of second 
order in displacement. 
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Basic Equations of Elastic Equilibrium in Vector Form 


The derivation of the equations herein is based upon the work presented by A. E. H. Love in 
Reference 6 and Landau and Lifshitz in Reference 5. Where possible, the notation used will follow 
that used by Love. 

Consider two adjacent cross-sections along the length of an elastically deformed beam that are 
a distance dz apart, as shown in Figure 25. Suppose also that the beam is deformed by forces that 
are applied along the length of the beam. The resultant effect that these forces have on the beam 
is estimated from: 

1, F their resultant force distribution measured per imit length along the longitudinal axis of 
the beam, and 

2. K their resultant couple distribution measured per unit length along the longitudinal axis of 
the beam. 

For the element shown in Figure 25 to be in elastic equilibrium, an internal stress distribution 
must be set up which will have the resultant effect of exactly counterbalancing the applied resultant 
force and couple distributions. 

Figure 26 isolates this element from the 
rest of the beam. Shown in this figure are: 

V = resultant of the internal stress 
distribution on the cross-section 
at 0^ 

V + dv = resultant of the internal stress 
distribution on the cross-section 
at o. 

These are the shear forces at o' and o respec- 
tively; similarly, 

iS = moment of the internal stress dis- 
tribution on the cross-section at 
o' and about the point o' , 

M + dM = moment of the internal stress dis- 
tribution on the cross-section at 
o and about the point o. 

These are the moments ato' and o respectively. 

For the element to be in equilibrium, the 
vector sum of the internal and applied force 




V +dV 


M + dM 


Figure 26— Beam element. 


31 


resultants must be equal to zero. Thus 


-V + V + dV + Fdz = 0, 


(67) 


or 


dz 


+ 


F 


0 , 


( 68 ) 


where dz is the distance between O' and 0 measured along the longitudinal axis of the beam. 

In addition to the above condition, the vector sum of the internal and applied couple or moment 
resultants must be equal to zero. Then 

- M + iii + diS + (- dzk) X (“ V) + icdz = 0 , (69) 


or 


m 

dz 


-;T-+kxV + K = 0 


(70) 


where 

k = unit vector normal to the cross-section at 0, and the moment due to the applied force 
resultant F about 0 is a second-order effect which can be neglected (see Love, Refer- 
ence 6). 

The vector equations that define the elastic equilibrium condition for any element along the 
beam are thus given by 


dV 

dz 


F 


0 , 


( 68 ) 


nM r - 

k X V + K = 0 . 


(70) 


For the equations to be effectively used in solving the problem of thermally-induced vibration, 
coordinate systems must be defined in which deflection can be measured and in which the applied 
forces associated with kinetic motion, damping, and the thermal environment can be defined. 
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Figure 27— Clamped free thin-walled cylinder of 
open section. 

Coordinate Systems 

Two coordinate systems are defined and are 
related to each other through an Euler angle trans- 
formation matrix. These systems are identical to 
those used in Reference 1. Consider the triads 
shown in Figures 27 and 28. 



Figure 28 — Boom cross-section at station z. 


Let; 

[Xj, Y^, Zj] be the inertially fixed triad defined at the centroid of the clamped root 

[x, Y, z] be the local, or body, triad. This triad slides with constant velocity along 
the centroidal axis of the boom, maintaining a fixed orientation with respect 
to the local cross-sectional principal axis of inertia (directions shown in Fig- 
ure 28) 

z be the coordinate used to measure arc length along the centroidal axis of the boom, 
measured positive from the clamped root. At z = 0, the inertial and local triads 
are coincident. 

[Tj, Tj, kj be a basis of three orthonormal vectors fixed in the inertial triad at z = 0 and 
parallel with the Xj, Y^^, and Zj , axes respectively, with their origin at the centroid 
of the cross-section 

[T, T, k] be a basis of three orthonormal vectors parallel respectively to theX, Y, and z 

triad at the point z, with their origin at the centroid of the cross-section at z = 0. 


The two sets of basis vectors are continually related to each other by means of an Euler angle 
transformation matrix. This transformation matrix is given below, and the rotation sequence is 
shown in Figure 29. 
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1 

0 

0 


cos 0^ 

0 

sin 0^ 


cos (p 

- sin cp 

— 

0 



0 

cos 0 j 

- sin ^ j 


0 

1 

0 


sin (p 

cos cp 

0 


J 

0 

L 

sin 0 ^ 

cos 0 ^ 

1 

- sin 0 ^ 

0 

cos <92 


0 

0 

1 


k 


cos 0 ^ cos cp 

sin 0 ^ sin ^2 ^^s cp + cos 0 ^ sin cp 
- cos 0 ^ sin0 ^ coscp + sin^j sincp 


- cos 0 2 sin cp 

sin<9j sin ^2 sincp+ cos 0 ^ cos cp 
cos 0 ^ sin 0 ^ sincp + sin 6 '^ cos :p 


sin 0^ 

- sin 0 ^ cos 0 2 
cos 0 ^ cos 0 ^ 


m 


( 71 ) 


L^J 


At any point z along the beam at time t, the unit vector which is tangent to the centroidal 
axis is defined by the vector equation 



Figure 29 — Euler angle rotation sequence. 


p, _ dX(z, t) 


(jY(z, t) 
dz 


(9Z(z, t) 


1<1.(72) 


where [x(z, t), y(z, t), Z(z, t)] are the coordi- 
nates of the point z along the centroidal axis at 
time t relative to the [Xj, Y^, zj axis system. 
Since k* is equal to k by definition and 


k - sin^^Cz, t ) i j - sin 0 j ( z, t)cos^ 2 (z. t)jj 

+ cos^^(2, t)cos^^2(^* > (73) 

it follows by equating components that actual de- 
flection can be obtained by solving the equations 


(9X(z, t) 
d z 


sin 6*2 (z, t) , 


(74) 


gY(z, t) 
d z 


-sin^j(z, t)cos 6 ' 2 (z, t), (75) 


^Z(z, t) 
d z 


cos j ( z , t ) cos 0 ^ ( z , t ) 


(76) 


Suppose that the body triad [X, Y, Z] slides along the centroidal axis of the boom with unit ve- 
locity. The angular velocity 0 with which the body triad rotates can be expressed in terms of the 
[i , T> k] set of basis vectors which maintains the same rotational orientation as [X, Y, z] relative 
to the inertial fixed-axis system. The angular velocity vector 0 can be written as 


n 


K\ + /< ' j + rk , 


(77) 
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where it caa be shown that 


d(9i d^2 

K - cos 0 ^ cos (p + sin (p 


d0. 




, - 1 ^ . ' 2 . 
k' - - COS ^2 Sin (p + cos (p , 


dz dz 


in . 


( 78 ) 

( 79 ) 

( 80 ) 


K and K * are simply the components of the curvature of the centroidal axis at the point z and 
time t, and r is the twist at the same point and time. The vector n will be needed in the subsequent 
analysis whenever the derivatives of i , j , and k with respect to the coordinate z are needed. 


Simplified Equations of Elastic Equilibrium 

The vector equations of elastic equilibrium as derived earlier are: 

dV 

^ + F = 0 , (68) 

^+kxv+i< = 0 . ( 70 ) 


In order to solve these equations, they must first be reduced to component form, be simplified by 
taking into account the constraints of the problem, and have the forcing functions defined. 

Using the notation and terminology of Love (Reference 6), we let: 

v = wT + N'i+Tk, ( 81 ) 

where N and N' are the "shearing forces," and the force T is the "tension"; 

M = Gi + G' T + Hk > ( 82 ) 

where G and G' are the "flexural couples," and the couple H is the "torsional couple." The forces 

N, N', T are called the "stress resultants," and the couples G, G', H the "stress couples." Further- 

more, the applied resultant force distribution is written as 

F = Fj T + F 2 T + F 3 k . (83) 
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and the applied resultant couple distribution is written as 


ic = Ki + K' j* + ©k ■ (84) 

By applying the principle that the change in a vector due to an infinitesimal rotation is equal 
to the vector product of the rotation vector and the vector itself, one is led to 







dz 


i = 

r j - /< ' k , 

(85) 

- 

dz 

0 X 

j = 

- ri + xk , 

(86) 





(87) 

dz 

n X 

k = 

^ ' i - /<j . 


The obvious substitutions can be made into Equations 68 and 70, the necessary vector operations 
performed, and the resulting expressions derived: 


fdN , , 1 


dN' 

1 


"8T ,, , , 

N' r+Tx' +fJ 

T + 

dz 

- T/< + Nt H-Fj 

T + 

+N'k+F3_ 


( 88 ) 


”80 

” 


dG' 


8H , , 

[dl- G' t + Hk' 

- N' + K 

i + 

-g^-H/<+GT+N+K' 

T + 

- G/<' +G' « + 0 


Small-displacement theory can be applied and the equations of equilibrium simplified if dis- 
placements from the unstrained state are assumed small. This assumption reduces the equations 
to a tractable form. The stability of the response of the system can now be analyzed; however, 
predicted steady state amplitudes of vibration must be interpreted with the assumptions of small- 
displacement theory in mind. 

By omission of terms that are second order in displacement, the equilibrium Equations 88 and 
89 become 


dN 


dN' ] 


'dT 

[dl + 

i + 

dz ^ 2 ] 

T + 

-J— F 
dz ^ 3 


(90) 


and 


dG 


fdG' ,1 


fdH 

[dF-N'+Kj 

i + 

l-dT +N + K-J 

T + 

[dT+ ®J 


k = 0 


(91) 
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From Timoshenko's theory of members of open section, the flexural and torsional couples are 
related to the components of curvature and twist by the following; 


G = EIk , 
G' = EI^<' , 


(92) 

(93) 


(9 2 r 


(94) 


where 

El = average bending stiffness, 

= torsional rigidity, 

= warping rigidity. 

If there is an initial pretwist in the rod defined by 

dcPp (z. 0) 
^0 ~ dz 


(94a) 


then the stress couple H is related to twist by 

H = - +C^(r-r^). (94b) 

It can be found in various texts, in particular Reference 1, that for a thin-walled cylinder of open 
cross-section 


where 


El = T [Ix + lj • 


Ct = 


hEr ^ ^ 


" . ^ 



sin ^ - 

r 2J 

1 cos -J 


(95) 

(96) 

(97) 


hr3 


1 . 

77 + ~2 j + *2 sini 


- (2v + 4>) hrh^ , 


(98) 
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(99) 


hr3 


7 + ^ • J 

l7T+ ■^1 - *2 sin0 


. 4> 

r sin "2 


77 


( 100 ) 


= shear modulus of elasticity, 

E = Young’s modulus of elasticity. 


Since small-displacement assumptions have been introduced into the equilibrium equations, 
they may be introduced also into the transformation matrix (71) to obtain 


i 
k 


COS cp sin cp 

“ sin cp cos qp 




-e^ 


~ ^ 2 ^ ^ 1 sin cp 

6 ^ sin cp + ^ j cos cp 
1 


r-. Y 


(101) 


into Equations 78, 79, and 80 (for curvature components and twist) to obtain 


d01 



dz 

cos ^ ^ 2 , sin cp , 

(102) 

d0, 

d0, 


6z 

sin cp + cos cp , 

(103) 


dcp 


r 

■ dz - 

(104) 


and into the displacement Equations (74, 75, and 76) to obtain 


. _ 5X(z, t) 

"^2 " c?z 

(105) 

N N 
>- 

1 

II 

Q3 

(106) 

N 

II 

N 

(107) 


The force equilibrium Equation (90) can be resolved into the inertial coordinate system, and 
the effect of tension deleted from the analysis. Equation 90 can then be rewritten as 



0 . (108) 
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Since the vector equation defining the sum of all moments (Equation 91) equals zero, each component 
must equal zero and hence 


dz 


dG' 

^ + N + K' = 0 

dH 


dz 


+ 0 □ 0 . 


(109) 

( 110 ) 
( 111 ) 


Substituting Equations 109 and 110 into Equation 108 yields 


d^G' dK' 


dz^ 


“dF + F J cos T 


d^ G dK 


dz 


2 dz '*'^2 


2 ^ sin cp 


d^ G^ 
dz ^ 


dK' 


s i n cp + ^ 


^d^G dK 

^ F^jcos cp 


0 , (112) 


By using the relation between G, G' and /<, k ‘ given by Equations 92 and 93, and by noting that 
products of small-displacement terms are of second order, we may express the derivatives of G 
and G' as 


d" G 
dz^ 


d^ G' 
dz ^ 


El 


El 


~ COS cp + 


sin ^0 


s i n rp + — CO s cp 


dz' 


dz^ 


(113) 

(114) 


Setting the components of Equation 112 equal to zero and substituting 113 and 114 yield 


d^ 

/dK' 


/dK 

dz^ 

\dz 



d^ 

fdjC 


, /dK 

dz^ 

y dz 

-fJ 



;in rp - 0 , 


(115) 

(116) 


From the relationships defined in Equations 94, 104, 105, 106, and 111, the equations of elastic 
equilibrium reduce to 

<? '* X( z . t ) /dK' 


"dz^ ^ td7 - FJ COSCP + 


(bT + Fj) sin 


(p = 0 , 


(117) 
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El 


(118) 


t) ^ 
d 


& " sincp 




cos cp 


0 , 


cp(z, t) t) 

dz^ T ^^2 


0 . 


(119) 


If pretwist is considered and if it and the amount of twist t about it can be assumed small- 
angle, then the torsion Equation (119) can be rewritten as 


[cp(z, t) - cpo (z, 0)] 

~ :~:r- + c, — 


[cp(z, t) - (Po (Z, 0)] 


+ 0 


0 . 


( 120 ) 


These are the equations of equilibrium which will be solved after the resultant force and couple 
distributions defined by F and ic have been derived. 


DISTURBANCE FORCES AND THE EQUATIONS OF THERMALLY-INDUCED VIBRATION 

The only disturbance forces which will be of interest in this paper are; 

1. The body forces associated with the translational and rotational acceleration of all mass 
elements along the length of the boom. 

2. The dissipative forces associated with viscoelastic damping of transverse and torsional 
motion. 

3. The thermal loading associated with the presence of the boom in direct sunlight. 

In addition to these forcing functions, the boom attached to an orbiting satellite will 
be disturbed by such effects as; solar pressure, aerodynamic drag, gravity gradient, 
satellite body motion, etc. When accounted for, each of these effects will perturb the 
actual motion slightly. By excluding these effects completely from the analysis, one can 
clearly show that thermally- induced vibration is possible and is a prime first-order effect 
not to be overlooked. 

The disturbance forces which are distributed along the length of the boom can be expressed 
as an equivalent resultant force and couple distribution. The equations of two-axis transverse 
bending and torsional vibration are the same as the derived equations of equilibrium under the 
action of the equivalent resultant force and couple distributions. 

The body forces associated with the reversed kinetic reactions of all accelerating mass ele- 
ments along the entire length must be obtained using vector techniques. The expression is not merely 
stated, but is derived from elementary considerations; since the standard smalL-angle assumptions 
break down when long, torsionally weak booms are studied. 
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The actual description of the dissipative forces is extremely complicated, and has not yet been 
completely defined experimentally. It is known (see Reference 10) that frictional contact in the 
overlap region is a strong contributor to the damping of large- amplitude transverse and torsion 
motion. It is also known that small-amplitude motion is damped by a viscoelastic-type damping 
mechanism. Experimental work done in a cold vacuum leads to results which give damping ratios 
for small-angle motion of anywhere from 1.0% to .01% of critical. No work, however, has been 
done to measure damping in a thermal environment or after extensive wear and/or thermal cycling. 

Since there are so many unknowns in the description of damping, it is assumed that by defining 
two viscous damping ratios, one associated with transverse motion and one associated with torsional 
motion, an estimation of the effect of transverse and torsional damping can be obtained. These two 
ratios are distinct from each other; no attempt is made to determine how or whether they are re- 
lated to each other through friction. They are simply used as independent parameters which can be 
varied independently of each other to determine their influence on the resultant motion of the boom. 

The thermal loading is dependent upon the boom's torsional response. By defining a sufficient 
number of thermal stations at which the torsional motion is recorded, and solving the thermodynamic 
equation at each station simultaneously with the vibration equations, a broken linear description of 
the temperature profile and hence the thermal loading is obtainable. By using higher-order inter- 
polation techniques between stations, an accurate description of the thermal loading can be derived 
in terms of a thermal bending moment and torque defined at every point along the length. 


Each of the three disturbance forces can be described completely in terms of a resultant couple 
distribution and a zero resultant force distribution. This is done in the following three subsections. 
In the fourth subsection, the results are combined and the actual equations of therm ally -induced 
vibration to be solved are written. 

Body Forces Associated With Kinetic Reaction of Mass Eiements 

The body forces associated with the reverse kinetic reactions can be derived from elementary 
considerations. At any particular instant of time, the vibrating boom will be in some deflected 
state, as shown in Figure 30. 

Let z = 0 define the clamped root of the boom, z = L define the free tip, z be an arbitrary 
point along the length, and z' be the location of any mass element between z and L. The mass ele- 
ment at z' is assumed to have mass p* dz' and rotational inertia l^dz', where 

p* = phP = mass per unit length of boom, 

I^=p*[r2 + e/ + 2e^sj = rotational inertia per unit length about the shear center axis of the 
boom. 

If the mass element at z' is separated from z by a finite amoimt and is accelerating in the 
transverse direction, it will produce a moment about the point z . 
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Ftgure 30— Insfanfaneous deflecfed position of boom. 


Let AR(z', z, t) define the vector connecting 
the points z and z ' . From the coordinate sys- 
tems defined, 

AR(z'.z,t) = [X(z' , t) -X(z, t)] Tj 

+ [y(z', t) -Y(z, t)] +(z' . (121) 

Furthermore, let p* R(z \ t) dz' be the force 
associated with the transversely accelerating 
mass element at z' ^ where 


(7 t “ at 


( 122 ) 


The moment at z due to the accelerating mass element is 


AR( z ' , z , t ) X p* R( z ' , t ) dz ' 


(123) 


and the equivalent couple distribution due to the reversed kinetic reactions associated with the trans- 
verse acceleration of all mass elements along the boom is where 


-r 


AR(z' , z, t) X p* R(z' , t) dz' . 


(124) 


Relative to the inertial coordinate system, K. can be written as 


dz 


-P* (2' - Z) 


Y(z' , t ) - 


+ p* (2' -2) 


d^X(z', t) - 


+ p" 


" j[X(z', t)-X(z, t)] - [y(z', t)-Y(z, t)] 


I a ^ x( z ' , t ) 


k, dz' , (125) 
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and by the transformation matrix (101) it can be expressed in the desired form relative to the local 
[i, j, ic] coordinate system as 



j Y(z', t) 

IL 


cos cp( Z , t ) 


X(z' , t) 
at2 


sin cp( z , t ) 


a^Y(z', t) 

<9t2 


sin cp( z , 


t) 


X(z' , t) 


COS (p(z , t ) 


- (z 


, • , 32 X(z, t) 52 Y(z', t) , ^ 32 Y(z, t)32 X(z', t) 

-z) — — + (z -z) 


dz^ 




dz^ 




k^dz',(126) 


where the differentiation has been carried out and second-order term deleted. 

Consider the kinetic reaction at z associated with the rotational acceleration of the mass ele- 
ment at z \ The torque at z due to this element is 


I 


s 


(9 ^ cp ( z ' , t ) 


k( z ' ) dz * , 


(127) 


where 

k( z ' ) = unit vector tangent to the boom at z ' . 

The equivalent couple distribution due to the reversed kinetic reaction associated with the torsional 
acceleration of all mass elements along the boom is where 


K, 



cp( z ' , t ) 


k( z ' ) dz ' . 


(128) 


To write relative to the body axes [i, j , k] , it must first be written relative to the inertial 
axis system , Tj > ] . Since 


k(z') - [k(z') ■ i J + [k(z') • j J Ti + • kj] kj , 

it can be seen from the transformation matrix (101) that 


= - I 


3 2 tp( z , t ) 
3t2 


3X(z, t) ^ 3Y(z, t) - ^ 

— 3i S + ~Tz — +'^1 


(129) 


(130) 


In terms of body axes, this becomes 


iCw 


cp(z, t) 



t) 


COS cp( z , 


t) + 


3Y(z, t) , , 

Jl sincp(z, 



+ 


dX{z, t) 
dz 


s Ln cp( z , t ) + 


(9Y(z, t) 
dz 


cos cp( z , 


t) 



• (131) 
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By combining Kg and K^ , the total couple distribution which describes the reversed kinetic re- 
actions of all accelerating mass elements can be written as K^,, where 


Kc = Kg + K^ 


K^ i + Kj,' j + 0J, k 


(132) 


and 


d^Y(z' , t ) 




COS cp( 2 , t ) 


d^X(z' , t) 


sin cp( z , t ) 


dz' 


t) 


ax(z, t) (3Y(z, t) 

cos<p(z, t) + sincp(z, t) 


, (133) 


r 


2 Y( z ' , t ) 5 2 X( z ' , t ) 

sincp(z, t) + cos rp( z . t ) 


dt^ 




dz ' 


<f(z, t) 
dt^ 


dX(z, t) (9Y(z, t) 

" sin(p(z, t) + ^ coscp(z, t) 


. (134) 




(z'-z) 


a ^ x( z , t ) Y(z' , t) _ a^ Y(z, t) a^ x(z ' , t) 

2 


dz^ 




dz^ 


dX.^ 


dz' - I 


rp( z , t ) 


(135) 


Dissipative Forces Associated With Viscoelastic Damping 

As previously discussed, experimental evidence shows that the dissipative forces are not pre- 
cisely definable for any arbitrary boom. It can be seen from Reference 10 that so-called identical 
booms may have damping characteristics which vary by as much as two orders of magnitude. In 
accordance with these findings, only viscous-type damping is included. It is argued that, by vary- 
ing the damping parameters, the effect of frictional damping can be estimated via an equivalent 
viscous damping term. 

Let R(z' , t) dz' be the dissipative force associated with the transverse motion, where 


/3 R(z',t)dz' = /3 


(9X(z', t) 
dt 


dz' 


dY(z* , t) 

+ /^x St dz- 


(136) 
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and let 


^(P(z\ t) 
^q> dt 


dz' 


(137) 


be the dissipative torque associated with the torsional motion. 

The equivalent couple distribution at z due to the dissipative force and torque associated with 
all moving mass elements along the boom is taken to be 


^ Z p 

^ = -^J (Z' -z)k, x,8^R(z'. 


„ a<p(z'. t) - 

0+/3cp kj 


dz' 


r 




^Y(z', t) - dX(z\ t) 


dt 




dt J 1 


(9cp(z, t) 




(138) 


By the transformation matrix (101), iCj, can be rewritten as 

K = i T + ©0 k . 


where, upon elimination of second-order terms, 




(9Y(z\ t) (?X(z', t) . 

cos (p(z, t) "ji sincp(z, t) 


dz ' 


f[^ 


^X(z', t) 

sincp(z, t) + ^ coscp(z, t) 


dz' 


0 


D 


drp(z, t) 
^cp dt 


(139) 


(140) 


(141) 


(142) 


Resultant Thermal Forces Associated With Thermal Loading 

In an earlier subsection, "Thermal Loading Associated with Non-Uniform Temperature Dis- 
tribution,^’ expressions are derived which define the resultant thermal bending moment and thermal 
torque at any particular thermal station z : If one assumes a sufficient number of thermal stations 
at which the thermodynamic equation is solved simultaneously with the vibration equations, then by 
interpolation techniques the resultant thermal loading can be accurately defined at every point along 
the length. 
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From the work previously done, the thermal loading at position z at time t can be defined by 
the thermal loading vector BM(z, t), where 


BM(Z, t) = BM^(Z. t) i + BMy (z. t) j + (z. t)k , 


(143) 


From Equations 64, 65, and 66, 


BM, (z.t) = 2 l^e^Ehr2p2 sin -2 2_^ 


n = 2 
n even 


(144) 


BMy(z,t) = -2i2e^ Ehr^ cos y 


qn(-. t) 


n - 1 
n odd 


( ri77r - P 


2 -p2 ’ 


Tsc (^- t) 


2t5e^EhrpJ^^ 

1 n odd 


_ ^ 
n77Pre^ cos p 

p2 - (ri'/7r)2 


d4^{ z , t ) 


t) d(p(z, t) 


^cp(z, t) dz 


(145) 


(146) 


To obtain the equivalent couple distribution 

Kh = K„r + k;T + 0„k 

due to the thermal loading. Equation 143 may be differentiated with respect to z to obtain the re- 
versed thermal reactions 


c?BM,, (z, t) 

Kh 


5BM (z, t) 

K„' = 


0 


H 


(^, t) 
dz 


where the derivatives of X, J, k are considered of second-order smallness. 
It must be noted that neither BM(z, t) nor 


dBM(z, t) 
dz 


46 



is necessarily equal to zero at the tip of the boom, z = L. This fact implies that the thermal load- 
ing introduces a time-dependent set of boundary conditions which must be satisfied for all time t . 


Thermoelastic Equations of Thermally-Induced Vibration 

Near the conclusion of an earlier subsection, ^'Simplified Equations of Elastic Equilibrium," 
it was found that the equations of equilibrium were 


El 


(9^ X(z, t) 
< 92 ^* 


-I- 


EI 


Y(z, t) 


+ 



= 0 , 


(117) 

(118) 


cp(z,t) a^cp(z,t) 

- + c, - . . ^ +0 = 0 ; 




az^ 


(119) 


and near the beginning of the same subsection it was found that the resultant force and couple dis- 
tributions F and K are given by 

j 

F = Fj i + F 2 T + F 3 k , (83) 


K = Ki + K' T + ©k . 


(84) 


For the three disturbance forces of interest, it was found that each could be described by a 
zero resultant force distribution and a particular resultant couple distribution. By vector addition, 
the total resultant couple and force distributions are 


F ° F^ i + F 2 j + F 3 k = 0 , 


(147) 


Kc + Kd + K„ 


(Kc + + K„) i + (K^' + V + K„')T ^ -(0c + ©D + ®„) k 


= Ki + K' T + 0k . (148) 

Upon substitution into the equations of equilibrium, it follows that 


El 


X(z, t) d , . . - d , . 

^ ^ +- 3 F (V +K„') COS9 +dF (Kc + K„ + K„) smcp = 0 


dz^ 


(149) 
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d ^Yjz, t _ ) ^ ^ ^ sin ^ - 4 - (Kc +K„ +K„) cos cp = 0 


(150) 


d^ q)(z, t) qp(z, t) 

C + C + 0+ 0+ 0 =0 

w (9z^ T ^z2 ’-'h ^ 


(151) 


where, if pretwist is considered, cp(z, t) in the torsion Equation (151) is replaced by [cp(z, t) 

-cPo (z. 0)] . By further substitution of the equations developed in the preceding three subsections 
and deletion of second-order terms, one obtains: 


a-*X(z, t) ^ ^ ^"x(z, t) ^ <^X(z, t) 

dz^ ^ at2 


BMy (z, t) BM^ (z, t) 

' “ coscp(z, t) + sincp(z, t) 


c?z" 


dz^ 


+ I. 


(9^ cp(z, t) Y(z, t) ^ (9^ cp(z, t) dY(z, t) 
^ (9z^ (9z (9t^ ^ 


(9f^ 


. (152) 


^ Y(z, t) . <92 Y(z, t) 

El ^^ 7 — - + /O* =^ 7 — - 

dz‘^ <9t2 


^Y(z. t) 

^ Px 9t 

d^BM(z, t) 


y. . , ^-'BM^Cz, t) 

— ^ sincp(z, t) — ^ coscp(z, t) 

^ dz^ 


dz 


<92 cp(z, t) (92 X(z, t) ^ (92 cp(z, t ) ^X (z, t) 


at- 


^z2 


t9 z (9 t ^ 5 z 


, (153) 


(9^ (P(Z, t) ^ q?(z, t) cp(z, t) t) 

dz' . (154) 

These are the equations of thermally- induced vibration which will be solved simultaneously 
with the thermodynamic equation defined at each thermal station. 

SOLUTION TO EQUATIONS OF THERMALLY-INDUCED VIBRATIONS 

Equations 152, 153, and 154 define the response of a thin-walled cylinder in a solar thermal 
field. These equations may have their solution expressed in terms of the orthogonal functions which 
define the undamped-unforced normal modes of vibration for the system. 


(z, t) 
dz 


f 


p” (z' -o 


92 X(z, t) (92 Y(z' , t) (92 Y(z, t) 92 X(z' , t) 


9z2 


9t^ 


dz^ 
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The boundary conditions that define the undamped-unforced system are time-independent. They 
define a boom that is clamped at the root, z = 0, in both two-axis bending and torsion. The tip, 
z = L, is free both to translate and to rotate, but has a tip weight attached to it which has both trans- 
lational and rotational inertia. The method of attachment of the tip weight is such that the tip re- 
mains free to warp. This condition is consistent with the end conditions forOGOIV,OGO V,andATSD. 

When the system is in a thermal field, there are certain time-dependent boundary conditions 
that must be satisfied and that are a function of the thermal loading. 

To solve a set of equations having time-dependent boundary conditions, the techniques de-? 
veloped by Mindlin and Goodman (Reference 11), and applied by Boley (Reference 2) to solve a dif- 
ferent thermally-induced vibration problem, are used. 

The resulting equations can be given a physical interpretation. They imply that the elastic- 
restoring force, which tends to drive the cylinder toward equilibrium, is measured not from its un- 
deflected position but from a time-varying position of static thermal equilibrium. The frequency 
content of the moving position of thermal equilibrium is therefore important, since if it has a strong 
component at the first natural bending frequency a resonant condition can be excited. 

To determine the motion of the position of thermal equilibrium, both the thermodynamic and 
the vibration equations must be solved simultaneously. 

Time-Dependent Boundary Conditions 

As previously discussed, the cylinder is assumed to be clamped in both two-axis bending and 
torsion at the root. Hence 


X(z, 

N 

11 

5 = 0 

- t)| 

(9X(z, t)| 

dY(z, t) 


dcp( z, t ) 

dz _ 

' z = 0 

dz 

z = 0 

dz 


(155a, b, c) 


(156a, b, c) 


In actual practice these root conditions are functions of the satellite body motion. For a gravity 
gradient satellite, the situation is fxirther complicated by the fact that roll, pitch, and yaw motion 
are coupled. Breakwell and Pringle* show that, for an undamped rigid body in circular orbit, low 
frequency (near orbital) non-linear resonances affecting gravity gradient stabilization can exist. 

This is another phenomenom not yet fully understood which can impair the performance of a flexible 
gravity gradient satellite. It is, however, not related to the relatively high frequency thermally in- 
duced oscillation problem discussed herein. 


At the tip warping is not restrained. Hence 


5 ^ cp( z , t ) 


z = L 


0 . 


(157c) 


If pretwist is to be considered, the boundary conditions on the torsion equation are obtained by 
simply replacing cp(z, t) by [cp(z, t) - (z, 0)] in the above and subsequent torsional boundary 

condition equations. 


♦Breakwell, J.V., and Pringle, R. Jr., ’’Nonlinear Resonance Affecting Gravity-Gradient Stability,” Astrodynamics (ed. Michal Lunc), 
Proc. XVI International Nautical Congress, Athens, 1965, Gauthier- Villair, Paris, pp. 305-325, 1966. 
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Attached to the tip is a tip weight which has both translational and rotational inertia. Let 
VI ^ = tip weight, lbs, 

= rotational inertia of tip weight, lbs sec^ in., 
g = acceleration of gravity, in./sec^ , 

and assume that the tip is free to translate and rotate in the thermal environment. At the tip there 
is zero bending moment attributable to the tip weight; but the thermal loading does impose a con- 
dition that must be satisfied. Hence 


" X(z. t) 
dz 


z - L 


j^BMy (Z, t ) cos cp(z, 


t ) + (z, t ) sin (p( 2 , t 


)] 


(157a) 


El 


d^Y(z, t ) 


Z - L 



t)sincp(z, t) - ( z'i t ) cos cp( z , 


■>] 


= L 


(15Yb) 


The translational and rotational acceleration of the tip weight gives rise to a shear force, and 
the thermal loading gives rise to what might be termed a "thermal shear." Both must be balanced 
at the tip by the internal stresses. Hence 


El 


X(z, t) 
dz^ 


X(z, t) _d_ 
1 g at2 C?Z 


j^BMy (z, t) cos cp(z, t ) 


+ BM 


(z, t) sin(p(z, t)j|- 


(158a) 


El 


Y(z, t) 


(9z^ 


) ^ 2 Y( z . t ) . d 


z. t) c/ r 

~ " -JI [BMy (z> t) Sin(p(z, t) 


BM 


^ (z, t ) cos (p(z, t)J 


(158b) 


+ c„ 


(9^ cp(z, t) c9cp(z, t) 


dz^ 


- C 


T dz 


^ (p(z , t ) 




(158c) 


These stated boundary conditions are consistent with the thermostress analysis of Boley and 
Weiner (Reference 2). Recent work by P. Jordan (Reference 12) shows that, when warping deforma- 
tion due to thermal strain is accounted for, the warping condition at the free end (Equation 157c) 
is non-zero. However, for a boom of practical interest (long length) this quantity is small and its 
effect approaches that of other deleted second-order terms. 
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Instantaneous Position of Thermal Equilibrium 


In order to solve the equations of thermally-induced vibration and satisfy the time-dependent 
boundary conditions, the concept of an instantaneous position of thermal equilibrium must be intro- 
duced. It is defined as the position that the boom would instantaneously assume if it had zero mass. 
From this definition, Equations 152, 153, 154, and the boundary conditions defined in the preceding 
subsection, it follows that the instantaneous thermal equilibrium shape of the boom is defined by the 
solution of the equations 

(Z, t) ^2 

^ t)coscp(z, t) + BM^ ( z , t ) s in q)( z , t )] , (159) 


El 


Y^(z, t) 
dz^ 



(z, t) sin(p(z, t)-BM^(z, 


t) cos cp(z, t)] , 


(160) 


(Z. t) <9%^ (z, t) <9T^^ (z, t) 


dz 




dz 


(161) 


and the boundary conditions 


(z, t) 

5Y^(z,t) 

5CPs (z, t) 

dz 

z = 0 ' '5Z 

-n ” 

z - 0 


= 0 


El 


a" (z, t) 


dz'^ 


El 


Y, (z, t) 


dz^ 


[BMy (z, t) cos cp(z, t) (z, t) sincp(z, t )] 

[BMy (z, t)sincp(z, t) (z. t) coscp(z, t)] 

(9^ cp (z, t) 


dz^ 


= 0 


EI- 


rz. t) 


dz^ 


El 


d^Y^(z, t) 


dz^ 


= [BMy (z, t ) cos rp(z, t ) + BM^ (z, t)sin(p(z, t)] 


= [BMy (z, t)sin(p(z, t) (z, t) cos cp( z , t)] 




dz3 


■“T dz 


= - t) 


(162a, b, c) 
(163a, b, c) 

(164a) 

(164b) 

(164c) 

(165a) 

L 

, (165b) 

L 

(165c) 
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Integrate both sides of Equations 159 and 160 twice with respect to z, integrate both sides of 
Equation 161 once with respect to z, apply the boundary conditions 164a through I65c, and find that: 


and 


El 


< 32 X 3 ( 2 , t) 


; - BM (2, t) cos (p(z, t) +BM^ (z, t)sincp(z, t) 

(9 z^ ^ 


El 


<32 Y3 (Z. t) 


BM (z,t)sin(p(z,t)“ BM^ ( z > ^ cp( z , t ) 


3^ 1>3 (2. t) dcf (2, t) 


<9z3 


'"t dz 


= - T3,(z, t) 


= Y3(z,t)|^.„ = cP3(z,t)|3.„ = 0 

5X3 (Z, t)| 


d2 


t) 

5Y3 (z, t) 

<3^3 (z, t) 


.=0 ' ^2 

z =0 ' ^2 


z = 0 


= 0 , 


<32 <P3 (Z, t) 


dz^ 


where it will be recalled from an earlier section that 


(166) 

(167) 

(168) 


(169a, b, c) 


(170a, b, c) 


(171c) 


4 > (2. t) , , 

BM^(z,t) = 2 Ehr^pz sin ^ 

n = 2 
n even 


BMj, (z, t) 


“ 2 V2 e Ehr^ cos 

c 


n = l 
n odd 


qn 

(n77r 


(65) 


T 


s c 


(z, t) 


2 i2e^ EhrP- 


CO 

L 

n = 1 
n odd 


n?7 


r. ^ 

n 77 Pre cos 

s ^ 

P^ - (oTrr 


P 

077- 


d^(z, t) 


t) dcp(z, t) 
(9{p(z, t) dz 


( 66 ) 
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Since the sim lineSL(z, t) is always assumed to be normal to the longitudinal axis of the boom, 
and (p(z, t) measures positive rotation of the cross-section about the axis, 

i//(z, t) = 0Q + cp(z, t) , (172) 


where 


'/’o - '/'(z. = o . 


Having defined a sufficient number of thermal stations at which the generalized thermal co- 
ordinate equations are solved, we have a good approximation to bm^ (z, t) , BMy (z, t), and (z. t) 
at hand. Equations 166 and 167 are initial-value problems and are readily integrated using any 
standard higher-order numerical integration technique. Equation 168 along with 169c, 170c, and 
171c defines a boundary value problem which is not as easily solved. It is, however, a linear dif- 
ferential equation with constant coefficients. 

To solve Equation 168, let fPp (z, t) be the particular solution of the initial-value problem de- 
fined by 


<5 ^ CPp ( Z , t ) 


dz^ 




t?cpp (z, t) 
dz 


(173) 


where 


T>P ( Z , t ) I 


t9(Pp (z, t) 
dz 


(9^ cpp (z, t) 
dz^ 


= 0 


(174) 


(this problem is readily solvable by numerical techniques). Making use of cpp (z, t) , it is then easily 
shown that the equation 


tPs (z. t) 


1 - cosh TjZ 


7 )^ cosh 77 L 


(Pp (L , t) 

dz^ 


+ tPp(z. t) 


(175) 


where 


V 


2 



salisfies Equation 168 and its respective boundary conditions. 


(176) 


53 



Thus the instantaneous thermal equilibrium position defined by the coordinates 


[x^ (z, t), (Z, t), (z, t)] 


can be obtained by numerically integrating the generalized thermal coordinate along the length of 
the boom at each instant of time. 

If pretwist is considered, cp^ (z, t) as it appears in the preceding equations is replaced by 
[cp^ (z, t) -(Pq (z, 0)] . The instantaneous position of thermal equilibrium is then defined by the 
coordinates [x^ (z, t), (z, t), cp^ (z, t) - (z, O)] . 

Equations of Undamped-Uncoupled Thermally-Induced Vibration 

The solution to the equations of thermally-induced vibration (152, 153, and 154) can be ex- 
pressed in terms of the solutions to the eigenvalue problem defined by the equations of undamped- 
uncoupled thermally- induced vibration; that is, 


El 


dz^ 


dt- 


dz^ dt^ 


dz^ 



[BMy (Z, t) cos qp(z, t ) + BM^ (z, t) sin cp(z, t)] , 

(177) 

t) ^ 

sin(p(z, t) - BM^ (z. t) cos cp(z, t)] , 

(178) 


cp(z, t) ^ ^2^(Z,t) _ Cz, t) 

(179) 

Ct 

5z2 = dt2 dz 

m by Equations 155a through 158c. Assume solutions of the form 

X(z, t) 

- X^ (z, t) + X(z) , 

(180) 

Y(z, t) 

“ y(z) . 

(181) 

cp(z, t) 

= (Z. t) + cp(z) , 

(182) 


and directly substitute into Equations 177, 178, 179 and the boundary condition Equations (155a 
through 158c), using the relations developed in the preceding subsection. One is immediately led 
to the following set of eigenvalue problems : 


El 


x(z) 

dz^ 


■ P* x(z) 


0 , 


(183) 
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El - P* “x^y(z) = 0 , 


(184) 


d‘* q>(z) 

dz'^ 


- C, 


cp( z ) 
dz^ 


" Is 'P(^) ■ 


(185) 


where 


x(z) 

l.=o = 

y(z)l z 

dx( z ) 


dy(z) 

dz 

z=0 

dz 

2 x(z) 


d^ y(z) 

dz^ 

z = L 

dz^ 

d^ x( 
dz3 


g 

y( 

dz3 

,'h. 

g 


^p(z)l 


_ dcp(z) 
dz 


0 , 


0 . 


d^ cp( z ) 
dz^ 


0 . 


0 ; 


cp(z) 

dz^ 


- C, 


dep(z) 
T dz 




0 . 


(186a, b, c) 
(187a, b, c) 

(188a, b, c) 

(189a) 

(189b) 

(189c) 


Equations 183 and 184 along with their appropriate boundary conditions are the standard eigen- 
value problem which defines the undamped-unforced normal bending modes of vibration of a clamped 
free beam with tip weight. Similarly Equation 185 along with its appropriate boundary conditions de- 
fines the undamped-unforced normal torsional modes of vibration of a thin-walled cylinder of open 
section without shear center effects, clamped at the root, free to warp at the tip, and having a tip 
inertia. The solutions to these equations are rather standard, but for completeness are included 
in Appendix B. 

From Appendix B the following is obtained. 

= nth natural frequency of vibration associated with transverse motion. 

n 

(z) = (z) = nth normal-'mode of transverse vibration; since bending stiffness El is con- 

stant about all axes, (z) and (z) are equal. 

Both (z) and (z) satisfy the orthonormality condition that 


f 


P" (Z)’'n (Z)'IZ - 




(190) 
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where 


P* (z) 



Wg = weight of boom without tip weight 

if m = n 

if m / n 

CO = nth natural frequency of vibration associated with torsional motion, 
cp^ ( z ) = n th normal mode of torsional vibration, 
cp^ (z) satisfies the orthonormahty condition that 

[ I, (^)'P„(^)<P„(^)dz = + 

*^0 



(191) 

(192) 

(193) 

(194) 


where 

0 < z < L 

(195) 

z = L 

If pretwist is considered, cp(z, t) and (p^ (z, t) are replaced by [cp(z, (z, o)] and [(p^ (z, t) 

-^Q (z, t)] respectively in Equations 179 and 182. 

Generalized Displacement Coordinates 

The results of the two immediately preceding subsections can be used to express the solution 
to the equations of thermally-induced vibration (152, 153, and 154) in terms of the normal modes of 
vibration and their corresponding generalized displacement coordinates. 

Assume a solution of the form 

00 

X(z, t) = X 3 (z, t) + 2^ (t)x„ (z) , (196) 

n = I 
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00 

Y(z, t) = Y, ( 2 , t) + b„ (t)y^ (z) 


(197) 


n = 1 


GO 

(P(z,t) = (P^(z,t)+ (t) q)^ (z) 


(198) 


Since at each instant of time the coordinates of the position of thermal equilibrium are known, 
they can be expressed as 


CD 


(199) 


Y, (z. t) 


00 


(t) y„ (z) 


( 200 ) 


n = 1 


( 201 ) 


where 


= Y. 


^ !>• 


(z)X^ (z, t) x„ (z) dz 


(202) 




f'- 

Jo 


(z) Y (z, t) y_^ (z) dz , 


(203) 


r„(t) 




I.C^ltPsCz. t)fPn (Z)dz . 


(204) 


Equations 196, 197, and 198 can now be rewritten as 


X(z 


CD 

■ t) = 2^ K (t)] ’‘n (^) 


(205) 


OP 

Y(z, t) = 7^ [;3„ (t) + b„(t)] y„(z) 


(206) 
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<p(z. t) 


(207) 


" 7^ [r„(t) + c„(t)]cp„(z) . 

1 


If pretwist is to be induced in the analysis, wherever cp(z, t) and cp^ (z, t) appear in the above 
equations the quantities [qj(z, t) (z, O)] and [cp^ (z, t) -cp^ (z, o)] should be substituted 
respectively. 

The equations of thermally-induced vibration (152, 153, and 154) are solved by: substituting 
into them the assumed solutions (205, 206, and 207); making use of the eigenvalue equations (183, 
184, and 185), which must be satisfied by all natural modes and frequencies of vibration; and mak- 
ing use of the equations (159, 160, and 161) that define the instantaneous positions of thermal 
equilibrium. Carrying through the prescribed algebra yields: 


n=l ^ 


+ p 


X (z) 


cp(z, t) Y(Z, t) ^ t) (9Y(z, t) 

dz^ ^ dzSt^ 


dt^ 


(208) 


r ^ K.(t)"b„(t)] +/3^^[/3„(t)+b„(t)] 


(t)^ (z) 


3^ <p(z, t) 3 ^ X( z , t ) ^ 3^ (p(z, t) 3X(z, t) 
3t^ 3 z 2 3z3t^ 


, (209) 


z: 


Is £-2 [yn(t) + C„(t)j +/3^^[r„(t)+C„(t)] 


Is % 'Pn (^) 


= !>• 


(z' - z) 


3^ X(z, t) 3^ Y(z* , t) 3 ^ Y( z . t ) 3X(z', t) 


3z^ 


3t = 


3z^ 3t^ 


dz' . (210) 


Equation 210 is unchanged if pretwist is included. 
Define the generalized displacement coordinates to be 
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where 


^ ( 211 ) 

Bn(t) = /3^(t) + b^(t), (212) 

c„(t) = r„(t) + C^(t). (213) 


Multiply both sides of Equation 208 by (z), Equation 209 by (z), and Equation 210 by cp^ (z) . 
Integrate over the length from z = 0 to z = l and apply the orthonormality relations given by Equa- 
tions 190 and 194 to obtain the ordinary differential equations which define the generalized displace- 
ment coordinates: 




dAp (t) 




dt ^ 

+ 2^ a> 

X 

n n 

dt 

+ K (o -°-p (t)] 

= Ap*(t). 

(214) 

B„ (t) 


dBp ( t ) 




dtz 

+ 2^ a; 

^ X X 

n n 

dt 

+ [Bp(t) -/9p (t)] 


(215) 



dCp (t) 




dt2 

+ 2^ 

cp 

dt 

+ [Cp(t)-7p (t)] 

= c;(t), 

(216) 


where A^* ( t) , b^* ( t) , and c^* (t) are. the generalized forces associated with the nth mode of vibra- 
tion, given by 




WT+W 3 


Is 


< 9 ^ cp(z, t) Y(z, t) cp(z, t) ,9Y(z, f) 

dz^ ^ azat^ 


Xp (z) dz 




Wt + Wb 


j 0 


Is (z) 


a^ cp(z, t) a ^ x( z , t ) ^ a^ cp(z, tj ax(z, t) 


at = 


dz‘ 


azat^ 5z 


y (z) dz 


(217) 


(218) 




= + 


It + Is 


I. I 


(z') (z' - z) 


~a ^ x( z , t ) a^ Y(z', t) _ a^ Y(z, t) a^ x(z',t) 


d z " 




(9 z^ 


(9t^ 


dz'dz .(219) 


and where the damping coefficients /5^ and^^ have been written in terms of their corresponding 
viscous damping ratios and If there is zero tip weight, they are defined by the expressions 

/3, = , (220) 

n 

/3<p = (221) 
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However, if there is a tip weight, since energy cannot be dissipated across it, its effect must be 
subtracted out. Thereupon it is found that 


r w 

(L)_ . 

^ ^ " "l^ + L (L) ■ 


( 222 ) 


(223) 


Generalized Forces 

Equations 217, 218, and 219 define the generalized forces A* (t), b^* (t), and (t) respec- 
tively. These quantities define the dynamic coupling between transverse and torsional motion. It 
is reasonable to assume that they can be adequately approximated by the first-mode bending and 
torsional components of the resultant motion. Accordingly, we let 


X(z, t) ^ Ai (t) xj (z) , (224) 

Y(z, t) = Bj (t)Yj (z) . (225) 

<p(z, t) = Cj (t) <pj (Z) , (226) 

and substitute into Equations 217, 218, and 219 to obtain 

d^Cj (t) 

V(t) = — (t)T , (227) 

C, (t) 

= -- ^ r- Ai(t)T„. (228) 


where 





d^ Yj (z) d(Pj (z) dYj (z) 

I. (^) 

-'o 1 

^1 dz dz 


X (z) dz 


(229) 


and 


c„*(t) = 


d^BjCt) d^Aj(t) 

— “-Bi(t) - 


dt^ 


(230) 
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where 


P 


n 


^^1’" 


(Z) 



(z') (z' -z) 


d2 Xj (z) 
dz^ 


( z ' ) dz ' dz . 


(231) 


It is noted that 


^n(^) - yn(^) 


(232) 


for all n. Since the equations for the mode shapes (z), (z), and cp„ (z) are known and given in 

Appendix B, the integrals defining and can be obtained by either closed-form or numerical 

integration. 

If numerical integration is used, then can be found by solving the differential equation 


dz 




I (z) 


^1 (z) 


d^ (z) dcpj (z) dy^ (z) 


dz 2 


dz 


dz 


X (z) 
n ^ ' 


(233) 


where 


'f'nCOL.o = 0 


(234) 


and 


Tn (^)l 


(235) 


Similarly, can be found by solving the simultaneous differential equations 

dr( 2 ) ^ 


dH(z) 

dz 


dT^ = H(z) 


“ /O* (z) y, (z) • 


dP„ (O 


dz 


1 d2xj(z) 

^ TT^TTl (-) r(z) . 


(236) 

(237) 

(238) 


where 


r(z)|^ = o - j* p(z) zyj (z)dz , (239) 
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z=0 


(240) 


H(2)| 



P(z) Ml (z) dz , 


P„(")l. = o = 0. 


(241) 


and 


Pn(^)iz = L " Pn- (242) 

It must be kept in mind that the generalized force expressions discussed in this section result 
from the fact that standard small-angle assumptions cannot always be applied in the analysis of ex- 
tremely long torsionally weak beams. For short and/or torsionally stiff beams, these forces should 
be negligible. 

By comparing the resultant responses obtained for various hypothetical examples with and 
without these terms included in the analysis, it has been found that: 

1. Results obtained either with or without A/ (t) and B/ (t) included are essentially the same 
for the class of problems being considered. For beams which are much weaker in bending, 
however, these terms could be important. 

2. Results obtained with and without C* (t) included are radically different for the class of 
problems being considered. As torsional rigidity is increased or length decreased, the 
effect of C* (t) on the resultant response is diminished. 

For production-t 3 /pe runs these observations are applied to increase the speed of digital com- 
putation. That is, c/ (t) , the coupling of bending into torsion, is left unchanged but all dynamic 
coupling of torsion into bending, defined by A/ (t) and B/ (t) , is assumed negligible and thus de- 
leted by setting 


= 0- (243) 

Summary and Method of Solution 

The work presented in the previous subsections shows that the solutions to the equations of 
thermally-induced vibration can be expressed in the form 


X(z, t) 


r. = 1 


A„ (t)x„ (z) , 


Y(z. t) = B„ (t)y„ (z) , 

n = 1 


(244) 


(245) 
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cp(z, t) 


(246) 


= 'Po(z.O) + C„(t)<P„(z)- 

n= 1 


The equations that define the generalized displacement coordinates are: 


d"A„(t) dA„(t) 

dt2 ~dt + K (t) (t)] 


d" Cl (t) 

• dt2 = 0 


(247) 


d^B„(t) dB^(t) 

+ 2C,_^ [B„ (t) (t)] = 


d^Ci (t) 


dt 


— Ai (t)T„ = 0 


(248) 


d" C„ (t) 
dt ^ 


dC„(t) 


+ 2Cq, % — [C„ (t) -r„ (t)] 


d"B, (t) d^AjCt) 

A.(o -irr--Bi(t)-^^ 


P„ • (249) 


Also shown is that the thermal loading defined by the quantities 

BM^ (z, t). BMy (z, t), (z. t) 


and hence the coordinates of the instantaneous position of thermal equilibrium defined by 


Xx (^. t) = 2^ “n (t)>‘n (") 


(199) 


Y3 (z. t) = 21! 

n = 1 


( 200 ) 




(201) 


can be written in terms of the generalized thermal coordinates (z, t) defined at a sufficient num- 
ber of thermal stations along the boom. 

To solve the generalized thermal coordinate equations 

(z. t) 

dt - Q„(Z, t), (250) 

the generalized thermal input (z, t) must be known as a fimction of time and position along the 
boom. 
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(z, t) has been shown to be an analytic function of the relative sun orientation ^(z, t). Thus, 
by simultaneously solving the generalized displacement and thermal coordinate equations, the gen- 
eralized thermal input Q„ (z. t) can be continually updated from a knowledge of 


■ACz. t) - ‘Aq + 9o 0 ) + (t)<p^ (z) , 

n=l 

and the quantities (t), (t), (t) can be updated from a continuous knowledge of (z, t). 

The numerical solution of the derived equations thus leads to the description of how a thin- 
walled cylinder of open section will respond when placed into direct simlight. 

The derived equations given here have been programmed for digital solution on the SDS-9300 
computer. The program could use as many as 20 thermal modes, 60 thermal stations, and 9 in- 
dependent modes of vibration. The results showed that good accuracy was possible using only four 
thermal modes, excluding that associated with . The number of thermal stations used for the runs 
excluding the effect of thermal torque was 20; when thermal torque was included, 50 thermal sta- 
tions were used. For all runs, transverse motion has been found to be essentially first-mode. Tor- 
sional motion, however, was found to have a slight second-mode component, so that two torsional 
modes of vibration were included in all runs to ensure a good representation of the thermal loading. 

Trapezoidal, Heun, Simpson^s, and fourth-order fixed-step Runge-Kutta integration techniques 
were used to solve the generalized displacement coordinate equation and in the integrations needed 
to define the thermal- equilibrium position. However, for production runs, Runge-Kutta integration 
proved the most efficient. To circumvent the numerical problems associated with simultaneously 
solving the thermodynamic and vibration equations when the thermal time constants are much 
shorter than the natural periods of vibration, an exponential interpolation technique was developed 
for solving the generalized thermal- coordinate equations. This permitted the use of a reasonable 
timewise integration step and prevented computational instabilities in solving the thermodynamic 
equations. 


PRESENTATION OF RESULTS 

The nominal model for the results presented is the antenna used on OGO IV and other existing 
satellites. These results should provide a heuristic understanding of the basic mechanism behind 
thermally-induced vibrations and the effect of changes in various parameters on its magnitude. The 
parameters studied are those which can be controlled by a manufacturer and used to govern the 
amplitude and stability of motion. Nominal magnitudes used to describe the geometrical and physi- 
cal properties of the antenna are shown in Table 3 . 

For all runs presented, the antenna is assumed to be initially in an unstrained state with a uni- 
form ambient-temperature distribution. The sun is assumed to be turned on at time zero, and its 
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Table 3 


Antenna material 

Antenna length 

Perimeter of cross-section 

Total overlap angle 

Wall thickness 

Weight density of material 

Shear modulus 

Young's modulus 

Bending stiffness 

Damping ratio for torsion 

Damping ratio for bending 

Torsional rigidity 

Warping rigidity 

Thermal conductivity 

Specific heat 

Solar radiation intensity 

Thermal expansion coefficient 

Emissivity 

Absorptivity 

Stefan- Boltzmann constant 


= beryllium copper, silver-plated 
= L = 60 ft 
= P = 2 in. 

= 0 = 90 deg 
= h = .002 in. 

~ P ~ .2714 Ib/in.^ 

= Gg = ,6 X 10^ Ih/ia? 

= E = .19 X 10^ lb/in.2 
= El = 16.6 lb ft2 
= ^ = .2 
= ^ .005 

= = .05 lb-in.2 

= = 1,700 Ib-in."^ 

= = 4.167 Btu/(hr-in.-°F) 

= c = .1 Btu/(lb-‘’F) 

= - 3.065 Btu/(hr-in.2) 

= = .104 X 10‘4 in./(in.-°F) 

= e = .035 
= = .13 

= cr = .121 xl0~^° Btu/(hr-in.^ ) 


orientation remains fixed in time, unless otherwise noted, along the inertially fixed negative axis. 
Throughout the ensuing dynamic response, the thermal input at all thermal stations is assumed to be 
that of one sun normal to the strained longitudinal axis. 

It should be further noted that in several examples the predicted amplitudes exceed the limits 
of small angle deflection theory. However, since numerical results agree with satellite data in a 
phenomological sense, the effects of including neglected terms would probably change the numerical 
values of the answers but not the basic characteristics of the motion. 

Effect of Thermal Torque 

Boley (Reference 2) shows that shear stresses exist in all cases where the temperature varies 
along the length of the boom. Appendix A derives the shear-stress distribution and resultant ther- 
mal torque expression theoretically for the special case of a thin-walled cylinder of open section. 
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Recent unpublished test data obtained by R. Predmore and C. Staugaitis at NASA/GSFC indicate 
that torsional response due to thermal torque can be reasonably well predicted by this analysis. 
More experimental work is needed, however, before a more positive statement can be justified. 

In addition to the work of Predmore and Staugaitis at NASA/GSFC, an independent experimental 
program has been carried out at NASA/ Ames by R, M. Beam. The results of Beam*s work have 
been published (Reference 13). Although Beam restricts his study to a boom with zero overlap, he 
has been the first to conclusively prove experimentally that the theory of thermally-induced vibra- 
tion as presented herein is physically realizable. In addition to demonstrating thermal flutter, 

Beam has in his report presented an extensive amount of experimental and theoretical work re- 
lated to the study of thermal torque. One of his most interesting experimental results is that, fpr 
certain sun orientations, thermal torque can excite a torsional instability. This region has been 
defined both experimentally and theoretically for a short zero-overlap boom. The material to be 
presented herein is consistent with the experimental findings both of Beam and of Predmore and 
Staugaitis. 

Keeping in mind the possibility of some inadequacy in any representation of thermal torque, we 
began a qualitative analysis of it by making numerous computer runs to determine the effect of 
changes in torsional damping, initial sun orientation, pretwist, length, and of the inclusion of the 
first three torsional modes of vibration. 

In order to build confidence in the analytic model, several of these computer runs attempted to 
duplicate laboratory test data. Actual data obtained by Predmore and Staugaitis at NASA/GSFC 
will not be presented herein. Simple visual observations concerning basic response characteristics 
will, however, be stated. 

For all booms tested, it was noted that, for some sun orientations relative to the outer seam of 
the boom, twisting was small; and that, for other sun orientations, twisting was quite large. To 
accent this effect, a 10-foot black painted standard boom was tested in the laboratory, and simu- 
lated with the analytic model. Because of the painting and end conditions, it was felt that damping 
could be reasonably well simulated by assuming 

" - 05 , = . 75 . 

Emissivity and absorptivity were measured and found to be 

e = .8 , = .8 . 

All other parameters were nominal. The effects of pretwist and relative sun orientation were stud- 
ied. This was done by slowly rotating the sun at 180°/min around the perimeter of the boom. The 
computed results are summarized in Figure 31. 

Tip twist <p(L, t) is plotted versus sun orientation 0(0, t) for varying amounts of pretwist. The 
curves labeled A, B, C, D, E, and F correspond to linear pretwists of 0, .5, 1.5, 3.0, 4.5, and 6.0 
degrees per foot, respectively. 
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Figure 3 1 -Thermally-induced fip fwist vs relative sun orientation for booms 

having linear pretwist. 

Both theory and test data show a significant change in the amount of tip twist when the sun 
moves across the outer seam of the boom. The magnitude of the twist observed runs approximately 
the same as that predicted for booms pretwisted 30 to 60 degrees. A significant dependence upon 
pretwist was not observed in the laboratory. It is suspected that, if shear center offset and the ex- 
act boundary conditions for the torsion equation discussed in Reference 12 were accounted for, the 
torsional response predicted would be less sensitive to changes in pretwist. 

In Reference 13 it is stated that if there is sufficient damping in the system all thermally- 
induced torsional motion will be stable. If damping is small, however, a torsional instability will 
exist for certain sun orientations. This phenomenon has been observed experimentally both by 
Predmore and Staugaitis at NASA/GSFC and by Beam at NASA/ Ames. Jordan (Reference 12) makes 
use of earlier work (Reference 1) by Frisch and predicts that the region of instability is where the 
function Ehr {dV/d^) shown in Figures 15 through 24 is negative. 

Beam, in his study of the zero -overlap boom, goes one step further and shows (by a small - 
angle single-degree-of-freedom analysis) that the instability can be either divergent oscillatory 
or divergent non-oscillatory. 

The results shown in Figures 32 through 37 attempt to illustrate the regions of stability and 
instability for a particular example. The example chosen is the black painted 10- foot boom used 
for the results presented in Figure 31. The transverse damping ratio is changed to 

= - 01 ' 

and a pretwist of 30° is assumed. The sun is turned on at time zero and starts to spin around the 
perimeter of the boom at a rate of ISOVmin in the positive direction starting from a position di- 
rectly over the axis of symmetry of the clamped root and over the overlap. 
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v>s(L,0 - TORSIONAL COMPONENT OF 

THERMAL EQUILIBRIUM POSITION «>(L,f) - TORSIONAL COMPONENT OF ACTUAL TIP POSITION 


t ‘ 


1 2 3 4 5 6 7 

TIME (min) 

Figure 32— Torsional component for tip of 10' boom, sun spinning at thermal torque included 




TIME (min) 

Figure 33— Torsional component of thermal equilibrium position for 10* boom, 
sun spinning at 180°/min, thermal torque included. 


In Figure 32 the amplitude of the torsional motion of the tip is plotted as a function of time 
Relative sun orientation and time are related by the simple linear equation 


0(0, t) - t , (radians) 


where t is time measured in seconds. 
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Figure 34— component for tip of 10' boom, sun spinning at 180®/m«n, thermal torque included. 



Figure 35— component of thermal equilibrium position for 10' boom, 
sun spinning at 180°/min, thermal torque included. 


Although the resultant response appears to be quite complex, it is evident that for certain sun 
orientations the amplitude of the torsional response grows rapidly and for other sun orientations 
it rapidly decays to an approximate steady state level. 

From this single figure it is not possible to pinpoint exactly regions of stability and instability; 
however, one does see that they exist and that they roughly correspond to those predicted by the 
small-angle analyses of Jordan (Reference 12) and Beam (Reference 13). That is, where the func- 
tion Ehr (dY/d^) shown in Figure 15, curve A,* is negative, linear theory applied to a single- 
degree-of-freedom model predicts a torsional instability. 

•180®/min - .0524 rad/sec ^ .05 rad/sec. 
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Figure 36— component for tip of 10* boom, sun spinning at thermal torque included. 



In Figure 33 the torsional component of the position of thermal equilibrium for the tip is plotted 
as a function of time. For the time span shown it is evident that after the initial thermal transient 
phase the curve shown repeats itself every two minutes (the time it takes the sun to complete one 
revolution about the perimeter). 

The points labeled (S) and (§) on this curve are the points at which a phase reversal occurs. That 
is, between these two points the motion of the torsional component of the thermal equilibrium posi- 
tion is such that it acts as a destabilizing negative spring. It has been remarked herein that 
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small-angle theory predicts an instability for relative sun orientations that yield negative 
Ehr (5v/(90) , A direct comparison between Figure 15 (curve A) and Figure 33 shows that the 
points (S) and (§) in the latter correspond to the points at which Ehr (dw/d^) equals zero in Fig- 
ure 15 and that for the time spans between @ and (§) the function is negative. 

As previously mentioned, the existence of torsional motion induces changes in the position of 
thermal equilibrium which in turn can induce transverse vibrational motion. This effect is shown in 
Figures 34 through 37, in which the and components of the actual position and thermal equi- 
librium position of the tip are plotted using the same horizontal scale as for Figure 32. 

It is evident that a very small motion of the position of thermal equilibrium at the first natural 
bending frequency of the system can and does induce a significantly strong response. A similar 
type of response has been observed experimentally. 

Although an easily grasped analogy between the phenomenon of thermal torque and some ob- 
vious everyday occurrence is not presently apparent, a heuristic discussion of this concept may 
prove useful for future investigation. 

Whenever a non-uniform temperature distribution exists along the boom, the normal thermal 
stresses that produce pure bending are accompanied by shearing stresses that produce torsion. 

For most closed- section tubes of practical interest, the shearing stresses are too small to pro- 
duce any noticeable twisting. For cylinders of open section, however, which are extremely weak 
in torsion, the effect of these shearing stresses can be quite significant. 

Since the mathematically perfect boom has never been made, non-zero shearing stresses will 
always exist in booms exposed to direct sunlight. The boom attempts to relieve this state of 
stress not only by bending, but also by twisting, toward a position of thermal equilibrium. In a 
single-degree-of-freedom analogy of torsion, the net effect of the shearing stresses (or their re- 
sultant thermal torque) is to act as a thermal spring whose magnitude and sign are dependent upon 
sun orientation. For the regions of torsional instability, the thermal spring is negative. 

It has been shown here that the inclusion of thermal torque in the analysis is necessary if a 
detailed description of the torsional response of the boom is desired. If one is interested only in 
whether or not transverse thermally-induced bending motion is stable, it has been found that the 
inclusion of thermal torque in the analysis of long booms (greater than 30') is merely a compli- 
cating factor that does not affect basic conclusions concerning the stability of transverse motion. 

Figures 38 and 39 are included to illustrate this point. The oblique (time) axis is introduced, 
and the origin of the coordinate system in which tip deflection is measured is allowed to simply 
translate along the time axis at a constant rate without rotation of the coordinate system. The mo- 
tion shown is only that of the actual position of the tip; neither the torsional response nor the posi- 
tion of thermal equilibrium toward which the boom is constantly moving is shown. 

In Figure 38 the effect of thermal torque is included in the analysis; in Figure 39 the effect is 
deleted. The boom studied is assumed to be 30' long and to have nominal physical and geometric 
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Figure 38— Tip mol-ion of 30' boom for nominal param- 
eters with thermal torque included when placed in di- 
rect sunlight at time zero. 


Figure 39— Tip motion of 30' boom for nominal param- 
eters without thermal torque included when placed in 
direct sunlight at time zero. 


characteristics. The sun is assumed to be fixed directly over the cross-sectional axis of sym- 
metry of the clamped root (along the negative Yj axis). In both figures the resultant response is 
seen to be unstable. Although the actual time histories do differ from each other, it can be seen 
that basic conclusions remain unchanged. 
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From this example and several other computer runs (not shown) it has been concluded that for 
long booms the effect of thermal torque is simply a complicating factor and that the prime inducer 
of torsional motion is the dynamic coupling of non-planar bending motion into torsion. 


By including only the dynamic coupling terms in the analysis, a practical answer to the stability 
question can be obtained. This is done in all succeeding runs, which show that with the effect of 


thermal torque deleted it is still possible to 
demonstrate that transverse thermally-induced 
vibrations can occur; and, furthermore, that the 
predicted amplitudes are reasonably close to 
those deducible from existing flight data. 

Effect of Length 

Figures 40, 41, and 42 illustrate the effect 
of a change in boom length by showing the com- 
puted dynamic response of a 60', 45', and 30' 
nominal boom respectively. For the time span 
shown, it is seen that each boom exhibits 
thermally-induced motion. The vibration ampli- 
tudes shown in Figures 40 and 42 are consistent 
with data obtained from the flights of the OGO 
IV and OGO V satellites: that is, the 60' boom 
on OGO IV is believed to oscillate with an am- 
plitude of about ±20 feet, and the 30' booms on 
OGO V are believed to oscillate with an ampli- 
tude of about ±2 feet. 

Besides showing this correlation between 
predicted and observed data, these three fig- 
ures also illustrate the sensitivity of the 
thermally- induced response to changes in an- 
tenna length. 

Effect of Transverse Damping 

Figures 43, 44, and 45 show the effect of a 
change in the transverse damping ratio of the 
antenna. If this parameter could be easily con- 
trolled and significantly increased above its 
nominal value of 0.005, the amplitude of the 
thermally-induced vibration could be stabilized, 
or at least limited to an acceptable magnitude. 



Figure 40— Tip moilon of 60* boom on OGO IV when 
placed in direct sunlight at time zero. 
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BOOM LENGTH 60 ff 
TRANSVERSE DAMPING RATIO 
0,005 


BOOM LENGTH 60 ff 
TRANSVERSE DAMPING RATIO 
0.05 


Figure 43— Tip motion of 60' boom for nominal param- Figure 44— Tip motion of 60* boom for nominal param' 

eters and transverse damping ratio .005 when placed in eters and transverse damping ratio .05 when placed ii 

direct sunlight at time zero. direct sunlight at time zero. 


Effect of Torsional Damping 

Figures 46, 47, and 48 illustrate the observation that, for large -amplitude thermally-induced 
transverse motion to exist, the antenna must respond torsionally at or near the first natural 


BOOM LENGTH 60 ft 
TRANSVERSE DAMPING RATIO 
0.10 


1 

10 ft 

Figure 45— -Tip motion of 60 ‘ boom for nominal param- 
eters and transverse damping ratio .10 when placed in 
direct sunlight at time zero. 


2.5 ft 

Figure 46— Tip motion of 30* boom for nominal param- 
eters and torsional damping ratio zero when placed in 
direct sunlight at time zero. 


transverse bending frequency co . This response in effect causes the position of thermal equilibrium 
to move at , and hence to drive the antenna at its resonant frequency. 

Equation 216 shows that, with the effect of thermal torque deleted by setting 7^ = 0, the tor- 

sional response depends upon the magnitude of the torsional damping ratio and upon the non-linear 
generalized force C* (t). The fact that (t) depends upon the actual position of the antenna and its 




Figure 47— Tip motion of 30' boom for nominal param- Figure 48— Tip motion of 30' boom for nominal param- 
eters and torsional damping ratio .10 when placed in eters and torsional damping ratio .20 when placed in 

direct sunlight at time zero. direct sunlight at time zero. 


position of thermal equilibrium leads to results that are extremely complex. The judicious choice 
of an example, however, can clearly illustrate the aforementioned observation. Accordingly, Fig- 
ures 46, 47, and 48 show how the resultant response of a nominal 30' antenna differs for the tor- 
sional damping ratios 0.0, 0.1, and 0.2 respectively. 
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In Figure 46, torsional damping is taken to be zero. From time zero, the torsional response 
shown in Figure 49 at the antenna tip gradually builds up at the first torsional frequency co^ . This 
has the dual effect of causing the position of thermal equilibriimi (shown in Figures 50 and 51) to 
change at about and of causing an averaging of the temperature distribution. Both effects are 
seen in Figures 50 and 51. Because the position of thermal equilibrium changes at the frequency 
, the effective thermal loading acts as a high-frequency forcing function. This has only a slight 



Figure 49—TorsionaI component of tip motion of the 30’ boom studied in Figure 46. 
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Figure 50— component of thermal equilibrium position of the 30' boom studied in Figure 46. 
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Y(L,t) -Yj COMPONENT OF X(L,t) - COMPONENT Y^(L,t)-Y, COMPONENT OF THERMAL 

ACTUAL TIP POSITION OF ACTUAL TIP POSITION ‘ EQUILIBRIUM POSITION 


effect on the transverse motion, as shown in Figures 52 and 53, and does not excite unstable trans- 
verse oscillation. The effect of temperature averaging appears in that the curves shown are closer 
to the time axis after transient response, rather than during it. 


In Figure 47, the torsional damping ratio is 0.1. Here the amplitude of the torsional re- 
sponse achieved during the initial thermal transients, shown in Figure 54, is seen to decay 



Figure 51— Y, componenf of the thermal equilibrium position of the 30' boom 



Figure 52— component of actual tip position of the 30' boom studied in Figure 46. 



Figure 53— component of actual tip position of the 30* boom studied in Figure 46. 
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gradually to zero at a non-resonant frequency. This decay causes the position of thermal equilib- 
rium, shown in Figures 55 and 56, to change at a non-resonant frequency and gradually decay to 
the position of static thermal equilibrium. This change causes the transverse motion, shown in Fig- 
ures 57 and 58, also to damp out to the displaced position of static thermal equilibrium. 
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Figure 54— Torsional component of Hp motion of the 30' boom 



Figure 55— component of thermal equilibrium position of the 30' boom studied in Figure 47. 



Figure 56— Y, component of thermal equilibrium position of the 30' boom studied in Figure 47, 
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Figure 57—X^ componenl- of actual tip 


position of the 30' boom studied in Figure 47. 



Figure 58— component of actual tip position of the 30' boom studied in Figure 47. 

In Figure 48, the torsional damping ratio ^ is 0.2. Here the torsional response, shown in Fig- 
ure 59, is found to increase gradually in amplitude, with a strong component at the first natural 
bending frequency . As the position of thermal equilibrium depends upon the torsional response, 
it too changes at about . This change is seen in Figures 60 and 61. The resultant effect is to 
drive the system at resonance and hence to excite the large-amplitude thermally-induced trans- 
verse motion shown in Figures 62 and 63. 


Effect of Torsional Rigidity 

Figures 64, 65 and 66 contain curves which indicate a practical means of designing a deploy- 
able antenna that will be stable in direct sunlight. The curves show the effect of varying the 
torsional rigidity of the antenna. Figure 64 shows the response of the nominal 60^ antenna, which 
is obviously unstable; in Figure 65, the torsional rigidity is increased by a factor of 10, to 
0.5 Ib-in.^. The slight increase in the response amplitude shows that, as the torsional response 
(not shown) has been reduced by a factor of nearly 10 in amplitude, the effect on tempera- 
ture averaging of large-angle torsional motion has also been reduced, with a corresponding increase 
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Figure 62— component of thermal equilibrium position of the 30' boom studied In Figure 48. 



Figure 63— component of actual tip position of the 30' boom studied In Figure 48. 


in the effective thermal loading. In Figure 66, the torsional rigidity is increased above its nominal 
value by a factor of 15, to 0.75 Ib-in.^ . Here the amplitude of the torsional motion has been re- 
duced to the point where it no longer has an effect upon the position of thermal equilibrium. The 
resulting transverse motion, as shown, is stable. 

Various manufacturers have used this idea and have submitted flight-quality antennas to NASA/ 
GSFC for evaluation; all the antennas submitted have the seams of the open- section cylinder zip- 
pered in various ways. Each of these designs exhibits a torsional rigidity at least 1,500 times 
greater than that of the nominal model defined, and about 100 times greater than needed to 
achieve stability. 
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BOOM LENGTH 60 ft 
TORSIONAL RIGIDITY I0 ft 
0.05 Ib-in2 


BOOM LENGTH 60 ft 
TORSIONAL RIGIDITY 


Figure 64— Tip motion of 60' boom for nominal param- Figure 65— Tip motion of 60' boom for nominal param 

eters and torsional rigidity .05 Ib-ln.^ when placed In eters and torsional rigidity .5 Ib-in.^ when placed 1 

direct sunlight at time zero. direct sunlight at time zero. 


Both RAE (Radio Astronomy Explorer), launched September 1968 and OGO-VI (Orbiting Geo- 
physical Observatory), launched June 1969, use the zippered- antenna concept. As predicted the 
four 750 foot long booms on RAE and the two 30 foot booms on OGO-VI are stable in and out of 
sunlight. 



Applications to Particular Spacecraft 

The original design specifications for sev- 
eral large spacecraft such as Nimbus, ATS, 
FR n, and AAP (Apollo Applications) have called 
for the use of long, tor sionally weak thin- walled 
cylinders of open section, some with and some 
without tip weights. The digital program used 
to derive the figures previously presented has 
been used to study the thermal response of each 
of these proposed boom configurations. The re- 
sults to be presented show that if these con- 
figurations were used a thermal instability 
could exist and seriously degrade the per- 
formance of the satellites. On the basis of 
these results, NASA has initiated a design 
change that will call for the use of zippered- 
type booms on each of the aforementioned 
satellites. 

The parameters (listed in Table 3, earlier) 
that define the boom used onOGOIV can be used 
to define the proposed booms on FR II, ATS, and 
AAP except for boom length and tip weight. 

For the French satellite, FR II, a 32’ boom 
had been proposed with a tip weight weighing 6 
lb and having a moment of inertia of .00258 
slug-ft^. The predicted response is shown in 
Figure 67. The predicted motion is unstable. 

An indication of the effect of introducing a 
tip weight into the system can be seen by com- 
paring this result with that shown in Figure 42. 
Both booms are of approximately the same 
length; however, one has a tip weight, and the 
other does not. It appears that the addition of 
a tip weight contributes to the instability of the 
system. 





Figure 66— Tip motion of 60 ‘ boom for nominal param- 
eters and torsional rigidity .75 Ib-in.^ when placed in 
direct sunlight at time zero. 


The proposed boom for ATS was 123.5’ long and had a tip weight weighing 8 lb with a moment 
of inertia of about .00266 slug ft^. The predicted response is shown in Figure 68. The predicted 
motion is unstable. 

The proposed boom for an AAP experiment was 120’ without tip weight. Again the predicted 
motion (shown in Figure 69) is unstable. The influence of tip weight on response characteristics 
can be seen by comparing Figures 68 and 69. 




Figure 67— Tip motion of 32' boom with 6-lb tip weight 
originally proposed for use on French satellite FR II. 



Figure 68— Tip motion of 123.5' boom with 8-lb tip 
weight originally proposed for use on ATS. 


The originally proposed boom for Nimbus differed significantly from that described in Table 3. 
Its deployment mechanism was designed so that the length of the boom could be changed by ground 
command. Preliminary studies indicated that lengths of 25', 35', and 45' would be used during 
flight; the stability of motion at each length was of interest. A list of the nominal parameters used 
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Table 4 


Uj 



Figure 69 — Tip molion of 120* boom with zero tip weight 
originally proposed for use on an Apollo Applications 
experiment. 




Antenna material 

Antenna length 

Perimeter of cross- 
section 

Total overlap angle 

Wall thickness 

Weight density of 
material 

Tip weight 

Tip inertia 

Bending stiffness 

Damping ratio for 
torsion 

Damping ratio for 
bending 

Torsional rigidity 

Warping rigidity 

Thermal conductivity 

Specific heat 

Thermal expansion 
coefficient 

Emissivity 

Absorptivity 


= polished beryllium copper 
= L = 45, 35, 25 ft 

= P = 3.5 in. 

= 0 = 155 deg 
= h = ,004 in. 

= P ^ .2714 Ib/in.^ 


= 10 lb 


- .002 slug ft^ 
El = 133.4 lb 


^(p 


= ^ = .005 


= .7162 lb-in.2 
- 3.546 X 10^ Ib-in.^ 


= 4.167 Btu/(hr-in.-“F) 

- c = .1 Btu/(lb-°F) 

- = .104X 10"^ in./(in.-°F) 
= e .06 

= a = .35 


to define the Nimbus boom is given in Table 4. The resultant response of the 25% 35% and 45^ 
booms is shown in Figures 70, 71, and 72. 


From these figures (70 through 72) it is evident that a boom which is significantly different in 
diameter, thickness, and absorptivity exhibits the same thermally -induced instability. It is recom- 
mended, therefore, that any thin- walled cylinder of open section be analyzed before it is used on an 
actual spacecraft. 
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CONCLUSIONS 

Most three-axis stabilized satellites de- 
ploying long booms have been plagued with 
varying degrees of anomalous spacecraft body 
motion. This motion, at times, is solely low 
frequency (i.e., near orbital), but in several 
instances, when the booms used are thin walled 
cylinders of open section, the motion is high 
frequency (i.e., near the fundamental beam fre- 
quency). With regard to the anomalous low fre- 
quency motion, there is a growing suspicion that 
the cause is a nonlinear resonance which allows 
exchange of energy between roll, pitch, and yaw 
motion. Since the boom studied in this report 
is, by definition, rigidly clamped in inertial 
space, the conclusions reached cannot be di- 
rectly applied to the low frequency phenomenon. 
With regard to the high frequency phenomenon, 
the following conclusions can be drawn; 

1. The theory of thermally -induced vibra- 
tion presented here can be used to explain the 
anomalous spacecraft-body motion observed on 
many three-axis-stabilized satellites deploying 
long, torsionally weak thin-walled cylinders of 
open section in direct sunlight. 

2. The problem of designing a thermally 
stable deployable antenna can be solved by zip- 
pering the seam of the open- section boom so as 
to increase its torsional rigidity by several 
orders of magnitude. 

3. Any satellite that requires long append- 
ages and geometric stability should not use tor- 
sionally weak thin-walled cylinders of open 
section. 
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Appendix A 


Thermal Bending Moment Components and Thermal Torque 


The derivations presented in this appendix are a brief summary of those given in more detail 
in Reference 1. 

The thermoelastic equation for longitudinal thermal stress at a point (s, z) at time t on the 
surface of the boom is 


(s, z, t) - e|€^ (s, z, t) “ [t(s, z. t) -T^]j , 


(Al) 


where 


(s. z, t) 

S (S’ 

e 

c 

T(s, z, t ) 

'V. 

T'a 

E 


longitudinal thermal stress, 
longitudinal thermal strain, 
thermal expansion coefficient, 
absolute temperature, 
absolute ambient temperature. 
Young’s modulus of elasticity. 


The longitudinal thermal strain of any fiber passing through an arbitrary cross-section along 
the boom can be adequately described by 


(S, Z, t) 


^0 (^’ (s) (Z’ t) + (s) (z, t) , 


(A2) 


where 


(s) = X- coordinate of point s on cross-section, 

(s) = Y- coordinate of point s on cross-section, 
iQ(z,t), fj(z, t),f2(z,t) = functions that are obtainable from the equilibrium conditions. 

The equilibrium conditions require that the resultant force and moments about the principal axes 
of the cross-section vanish; that is. 


hcr^ ( s , z , t ) ds 

•^0 


ho-j (s, z, t) (s) ds 


hc7^ (s,z, t)x^(s)ds - 0 


(A3) 
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It follows from the determination of the three functions f^Cz, fj(z, t), and (z, t) that 
the longitudinal thermal stress is 


BMy (z, t) BM^ (z, t) 

(s,z. t) = Ee^ [t^ (z, t) -T(s, z, t)] j (s) + — (s) 


(A4) 


where 


= r 


hy^^ (s) ds 


(A5) 


•'o 


ly - 1 hx/(s)ds. 


(A6) 


-'c 


■ P t) ds , 


(A7) 


r 


BiM^ (z, t) - Eh T(s , z, t ) ( s ) ds 


(A8) 


BMy (z, t) 


•'o 


T(s, z, t) X (s) ds. 


(A9) 


The quantities bm^ (z, t) and bm^ (z, t) are the thermal bending moment components about the x 
and Y principal axes of inertia, respectively, of the cross-section at(z, t) . 


For thin-walled members of open section, the shear stress t(s, z, t) is related to the longi- 
tudinal stress cr^ (s, z. t) by the partial differential equation 


c)hT(s, z, t) _ , '•) 

ds ^ ^ 7z 


(AlO) 


Substituting Equation A4 along with the geometric identities 




y<;(S) - - 



(All) 


(A12) 
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into Equation AlO, applying the chain rule of differential calculus to eliminate differentiating the 
temperature distribution with respect to z, and integrating both sides of the equation from 0 to s 
yield the expression that defines the shear stress distribution. That is, 


hr(s, z, t) 


L 


r 2 <9BMy (z, t) p 

+ t;; d^p{z. t) ■ \T ■ Tj T 


r 2 <9BM^ (z. t) 

t; a^(z, t) 



0 

sin ~2 



d\p{ z, t) dq)(z, t) 
‘(9cp(z, t) ~Sz 


(A13) 


The element of torque (z. t) due to the shear stress at an element ds of unit longitudinal 
length on the surface can be obtained from the expression 

dTg^ (z,t) = “^sc^ hr(s, z, t) ds , (A14) 


where 

= vector from the shear center to the surface element ds. 
Since may be written as 


(A15) 


where 


= vector from the shear center to the geocenter of the cross-section, 
T - radius vector to the surface element ds, 

it follows from geometry that 


|e^ xhr(s, z, t)| 


e ^ hr( s , z , 


t ) cos 




(A16) 


I F X hr( s , z , t ) I - rhr( s , z , t ) . 


(A17) 
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Hence the thermal torque is obtainable from the expression 


T (z, t) = - e 

*sc V ’ / s 


r /s 4 >\ 

j hr( s , z , t ) cos I hr( s , z , t ) ds , 


(A18) 


where hr(s, z, t) is defined by Equation A13. After the appropriate substitutions and integrations 
have been made, the expression for the thermal torque is found to be 


T 


sc 


(z. t) 


dV( z , t ) t) dcp( z , t ) 

^ (9i/^(z, t) ^cp(z, t) dz 


(A19) 


where 


V(z. t) 


^+-1 r 

2^2 

J Jo 


T(s , z , t ) ds 



r^s /s 1 

- - 2 ) + 1 


0 

\ / 



T(^, z , t ) d^ ds 


(A20) 


and 



(A21) 
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Appendix B 


Uncoupled Bending and Torsional Modes of Vibration 


In the subsection, ^Equations of Undamped- Uncoupled Thermally -Induced Vibration," it is 
shown that, to obtain the solution to the equations of thermally-induced vibration, the uncoupled 
normal modes of bending and torsional vibration must be determined. These modes of vibration are 
the linearly independent solutions to the eigenvalue problems defined by the differential equations 


d"* x(z) 

El P x(z) - 0 


(Bl) 


El 


y(z) 

dz'^ 


p* y(z) 


(B2) 


d"^ cp( z ) d^ cp( z ) 


dz'^ 


■"T 




(B3) 


and their respective boundary conditions 


X(0) = y(0) ^ cp(0) 


dx(0) _ dy(0) d(p(0) 


dz 


dz 


dz 


(B4a, b, c) 
(B5a, b, c) 


x(L) _ d^ y(L ) ^ d^ cp(L) ^ 
dz^ dz^ dz^ 


(B6a, b, c) 


d^ x(L) 
EI-^ 


= 0 


(B7a) 


El 


d^ y(L) 

dz^ 


y(L) = 0 


(B7b) 


d^ cp(L) dcp(L) 


dz^ 


Ct - It^cp = 0 . 


(B7c) 
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Since Equations Bl and B2 are identical in form and have the same boundary conditions, only the 
solutions to Bl and B3 will be outlined. 


Bending Modes of Vibration 
Let 




p* a; 2 

^ X 

“ El“ 


(B 8 ) 


Then Equation Bl may be rewritten as 


x( z ) 


- x( z ) 


0 . 


(B9) 


The general solution to this equation is 


x(z) - Dj cosh \z + D2 sinh\z + D3 cos Xz + sin\z . 


(BIO) 


Directly substituting into the four boundary -condition Equations B4a, B5a, B 6 a, and B7a yields 


0 = + D 3 , 


0 = Dj + , 


0 = Dj cosh \L + Dj sinh\L “ Dj cos X.L " sin\L 


(BID 

(B12) 

(B13) 


'^T \ / '^T 

0 - ^sinhX.L + XL cosh \Lj + D2 I cosh XL +'^XLsinhXL 




+ D, ( s in XL + XL cos XL ) + I " cos XL + XL s in XL 






(B14) 


where 

W 3 = p*Lg . (B15) 

For a non-trivial solution of these equations to exist, the determinant of the coefficients of Dj , , 

D 3 , and must be equal to zero. The evaluation of this determinant yields what is commonly 
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referred to as the "frequency equation." It is found to be 


_ 1 -r sec KL sech \L 

XL (tan XL~ tanh \l) 


(B16) 


The nth solution to this equation (B16) in ascending order of magnitude will be written as 

It follows from Equation B8 that the nth natural frequency of the clamped free beam with tip weight 

is 



(B17) 


Associated with each eigenvalue (^L)^ there exists an eigenvector. This eigenvector is the 
nth mode of vibration. It is obtained by substituting 


X 

n 


L 


(B18) 


into BIO to obtain 




D, cosh X 

In n 


sinhX 

2 n n 


Z + 


D_„ cos X 

3 n n 


Z -i- 


s in X z . 

4n n 


(B19) 


Directly substituting B19 into the boundary- condition Equations B4a, B5a, and B6a and solving the 
equations simultaneously yield 


X ( z ) = D, cosh X z 

n In n 


sinX z sech X z “ tanh X z 


s in X^ L sech L + tanh X^ L\ 

cLx LsechX L + 1 y (cos X„ z sech X„ z - l) 


(B20) 


In order to define uniquely the magnitude of Dj„ , the modes must be normalized. The ortho- 
gonality relation that defines Dj„ for this analysis will be taken to be 


f 


P’ (Z)dz 




(B21) 


where 





0 < z < L 


z = L 


(B22) 

9P 



0 m r n 


To7'sio}?al Modes of ViM^ation 

For numerical computation and notational purposes it is convenient to solve Equation B3 in 
terms of dimensionless quantities. 


then Equation B3 becomes 


9(2* ) cp(z* ) 


— c ~ — ) 


and the boundary conditions B4c, B5c, B6c, and B7c become 


_ dcp(O) _ d^(p(l) _ . 


(B27, B28, B29) 


d3<p(l) d^,(l) "(P 


Equation B26 has the general solution 


cp(z*) - Ej sinhkaz* + coshkaz* + sinkyz* -s- E^ coskyz* , 


where 


I. 
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r = 



(B34) 


(B35) 


In a manner analogous to that of the preceding subsection Equation B31 is substituted into the four 
boundary- condition equations B27 through B30 to yield 

0 = Ej + (B36) 

0 = kaE, + kyEj . (B37) 

0 - Ej sinh ka -r cosh ka - k“ sin ky ” E^ k^ y^ cos ky , (B38) 


0 - cosh ka + j ^ k ^ s inh ko^ 




^ ^2 (” sinh ka -r ^ \^2 ^osh ka I 


+ ( k^ ya^ cos ky -i- ^ X.^ sin ky 


E^ ( - k^ ya^ s in ky cos ky ) . (B3 9) 


)• 


For a non-trivial solution of these equations to exist, the determinant of the coefficients of Ej , £3 , 
E 3 , and must be equal to zero. The evaluation of this determinant yields the frequency equation 



2X.^ sec ky sech ka + 



+ k\ t an ky t anh ka 


\k2 


t SP" 


(a tan ky -y tanh ka) 


(B40) 


If the substitutions 


X *2 = \2 ,^2 


(B41) 
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are made, the frequency Equation (B40) may be rewritten as 


I L 


2 



sech ka* sec ky* 



+ /V* tan k7* tanhka* 


kK* 





(a* tanky* -y* tanhka"') 


(B42) 


(B43) 


(B44) 


It is obvious that either Equation B40 or B44 can be used to determine the natural frequencies 
of the cylinder of open section with tip inertia when the tip is free to warp. It is convenient when 
numerically solving the frequency equation, however, to use Equation B40 for cases where k^ > 100 
and to use Equation B44 when < 100. This choice is made to avoid the inherent problems associ- 
ated with digitally taking the difference between very large nearly equal numbers. 

Let be the nth value of k that satisfies the frequency equation in ascending order of mag- 
nitude. It follows then that the nth torsional natural frequency as defined in Equation B32 is 

C \ 2 

0.2 = (B45) 


The nth eigenvector, that is, the nth torsional mode of vibration, is obtained by substituting into 
the general solution given by Equation B31. This yields 




E, sinhko. z* cosh ka z* 

In n 2 n n 


E^ sinky z* 

3 n ' II 


E^ COS ky 2 * , 

4n 'n 


(B46) 


where 


I, 


\2 


a 


n 


7 


n 



(B47) 


(B48) 


(B49) 
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Substituting into Equations B27, B28, and B29 and solving for the unknown coefficients yield the nth 
torsional mode of vibration. This is 


<P„ (z) - E, 


ka^ 

cosh — z 


ka^ 


tanh 


z “ sech 


ka 

n 

L 


w 


where 


/ \" 
G„ - sech — z cos — zj 


(B50) 


= z* L , 


a (a tanh ka +y sech ka sin ky \ 

_ n\n n'n n ^ I 

a 2 + y 2 ka cos ky 


(B51) 

(B52) 


The magnitude of E,„ can be specified by requiring that the mode shape satisfy the normaliza- 
tion condition 


f I. (zlVn 

Jo 


(B53) 


where 


I. (z) 



0 < z < L 
z = h 


(B54) 
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Appendix C 

List of Symbols 


Lowercase English 


( t) amplitude coordinate associated with the nth transverse bending mode in the 

(Xj, Zj) plane 

(t) amplitude coordinate associated with the nth transverse bending mode in the 

(Yj.Z^) plane 

(t) amplitude coordinate associated with the nth torsion mode 


c 


specific heat of material 


thermal expansion coefficient of material 

distance between the geometrical center and the shear center of the cross- 
section 


g 

h 

[i, T. k] 


[ij. Ti. kj] 


m 


acceleration of gravity 
thickness of cylinder wall 

orthonormal set of basis vectors parallel to body axes triad [x, Y, z] re- 
spectively at (z, t), but having their origin at z = 0 

orthonormal set of basis vectors parallel to inertial axes triad [x^ , Yj , z J 
respectively, having their origin at z =0 

positive integer 

positive integer 


n 


positive integer 


n(s) 


unit vector normal to the cross-section and directed outward at point s on 
the circumference 


(z, t) nth generalized thermal coordinate at thermal station z and time t 

r radius of cross-section 
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s 


arc length measured from outer to inner seam around the cross-section 


t coordinate used to measure time 

(z) nth normalized transverse bending mode of vibration in the (x^ , plane 

(s) distance from the point s on the cross-section to the Y-axis 

(z) nth normalized transverse bending mode of vibration in the (y^ , ) plane 

(s) distance from the point s on the cross-section to the X-axis 

z coordinate used to measure arc length along the centroidal axis of the boom 

z ' arbitrary point on centroidal axis between z and L 


Uppercase English 


A„(t) 

A„* (t) 

B„(t) 

B;(t) 

BM^ (z, t) 
BMy (z, t) 
Cn(t) 


Ct 

Cw 

Di 

E 


generalized displacement coordinate associated with the nth transverse bend- 
ing mode in the Xj , Zj plane 

generalized force associated with the generalized displacement coordinate 

An(^) 

generalized displacement coordinate associated with the nth transverse bend- 
ing mode in the (Y^ , z^) plane 

generalized force associated with the generalized displacement coordinate 

B„(t) 

thermal bending moment about the X-axis at (z, t) 
thermal bending moment about the Y-axis at (z, t) 

generalized displacement coordinate associated with the nth torsional mode 

generalized force associated with the generalized displacement coordinate 
Cn(t) 

torsional rigidity 
warping rigidity 

constants of integration (i = 1, 2, . . . ) 

Young’s modulus of elasticity 


El 


bending stiffness 
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F 

[F,. F2.F3] 


G 

G' 


Gs 

H 

I 


Js 

K, K' 


Kb’ K3' 

Kc, K,' 

Kb’ K3' 
Kh« K; 
Kw k; 

Kt 

L 

M 

N 

N* 


resultant force distribution along boom 

components of resultant force distribution with respect to basis vectors 

[i , T> k] 

component of stress couple parallel to the x-axis at (2, t); ’’flexural couple” 
component of stress couple parallel to the Y-axis at (z, t); ’’flexural couple” 
shear modulus of elasticity 

component of stress couple parallel to the z-axis at (z, t); ’’torsional couple” 

average cross-sectional geometrical moment of inertia 

cross-sectional geometrical moment of inertia about the x-axis 

cross-sectional geometrical moment of inertia about the Y-axis 

mass rotational moment of inertia per unit length about the shear center axis 

mass rotational moment of inertia of the tip weight 

solar radiation intensity 

components of resultant couple distribution per unit length parallel to the 
X-, Y-axes respectively at (z, t) 

components of K, K' associated with kinetic reaction of transversely acceler- 
ating mass elements, respectively 

sum of Kg, and respectively 

components of K, K' associated with the dissipative forces 

components of K, K' associated with the thermal loading 

components of K, K' associated with kinetic reaction of torsionally accelera- 
ting mass elements 

thermal conductivity of material 

antenna length 

moment of the internal stress distribution in the cross-section about its 
centroid 

component of stress resultant parallel to X-axis at (z, t) 
component of stress resultant parallel to Y-axis at (z, t) 
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p 

Q„(z. t) 
R(z, t) 

^(z, t) 

T 

To 

T(s, z. t) 
T(s, z, t) 
Tn(0 
Tsc("> 

V 

V(z, t) 

[X, Y, Z] 

X 


total perimeter of cross-section, tape width 

generalized thermal input associated with the nth thermal mode and nth gen- 
eralized thermal coordinate q„ (z, t) at (z, t) 

transverse acceleration of the mass element at z at time t 

unit vector normal to z-axis and directed toward sun at (z, t) 

component of stress resultant parallel to Z-axis at (z, t); "tension^' 

steady state mean temperature of cross-section in thermal field 

absolute temperature of the point (s, z) on the surface of the boom at time t 

deviation of the absolute temperature T(s, z, t) from the mean 

nth normalized thermal mode 

thermal torque at (z, t) about axis normal to cross-section 

resultant of the internal stress distribution on the cross-section 

thermal torque coefficient used to define (z, t) 

boom weight 

tip weight 

body axes at (z, t) 

body axis at (z, t) ; principal axis of inertia normal to axis of symmetry and 
longitudinal axis, directed so that x, Y, Z is a right-handed system 


Y body axis at (z, t); principal axis of inertia parallel to axis of symmetry 

and normal to longitudinal axis, directed away from overlap 

Z body axis parallel to longitudinal axis, directed positive from the point z 

toward the tip of the boom 

[Xj , Yj , zj inertially fixed axis system, coincident with the body axes at z = 0 

[x(z, t), Y(z, t), z] displacement coordinates of the point z along the centroidal axis at time t in 
the [x^, Y^, zj coordinate system 

[x^ (z, t), Y^ (z, t)] displacement coordinates of the point z along the centroidal axis of the boom 
when it is aligned along the instantaneous thermal equilibrium shape at time 
t measured in the [Xj , Y^ , Z^ ] coordinate system 
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Script Greek 




^„(t) 

^<p 

>„(t) 

s 

c 

5 

m, n 

e 

^x 

^x 

n 

^q) 

V 


p 

p* 


(J 


T 


T 

n 

CP(2, t) 


amplitude of nth mode at time t in the modal expansion of (^, t) 
absorptivity of boom surface 

amplitude of nth mode at time t in the modal expansion of (z, t) 

viscous damping coefficient for transverse motion 

viscous damping coefficient for torsional motion 

amplitude of nth mode at time t in the modal expansion of (z, s) 

distance between geometrical center and centroid of cross-section 

Kronecker delta function 

emissivity 

transverse damping ratio 

transverse damping ratio of nth transverse bending mode 
torsional damping ratio 

torsional damping ratio of nth torsional mode 
defined constant = yc^/C^ 

Euler angle sequence used to define rotational orientation of [ i , j , i<] rela- 
tive to [ ij , Tj . kj 

curvature component about the x-, Y-axes at (z, t) 

thermal decay constant of nth thermal mode 

weight density of material 

mass per unit length of boom 

Stefan- Boltzmann constant 

twist about z-axis at (z, t) 

thermal time constant of nth thermal mode 

angular amoimt of overlap 

angular amount of twist at (z, t) of the cross-section about its longitudinal 
axis 
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(2’ t) 


i/f(Z, t) 


torsional coordinate of the point z along the boom when it is aligned along the 
instantaneous thermal equilibrium shape at time t 

nth normalized torsional mode of vibration 

angular position of the sim measured relative to the cross-sectional principal 
axes of inertia 

nth natural frequency of transverse vibration 
nth natural frequency of torsional vibration 


Printed Greek 


B 

r 

z, t) 
E 
H 
0 


ic 

E(s, z, t) 


T 


n 


functions used in defining generalized thermal input (z, t) 
defined parameter 

vector connecting points z‘ and z along the centroidal axis at time t 
functions used in defining generalized thermal input (z, t) 
defined parameter 

component of resultant couple distribution per unit length parallel to Z-axis 
at (z, t) 

component of 0 associated with kinetic reaction of torsionally accelerating 
mass element 

component of 0 associated with the dissipative force 

component of 0 associated with the thermal loading 

applied resultant couple distribution along boom 

defined to distinguish between sunlit and dark sides of boom 

parameter used in generalized force expression; associated with nth tor- 
sional mode 

parameter used in generalized force expression; associated with nth trans- 
verse mode 


T(s, z, t) function used to define the heat input about the surface point (s, z) at time t 

nth eigenvalue of the thermal characteristic equation 

n curvature vector; angular velocity vector expressed relative to [T , T> k] of 

the [X, Y, Z] axes when it slides at a constant rate along the longitudinal axis 
of the boom 
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